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_ Abstract. The paper gives a solution ‘of the problem of motion of a particle in the potential field 
V = —1/r + 2kP2 (sin 0)/r? + k’Ps (sin 6)/r5, 


where k and k’ are small parameters. The first approximation is furnished by a Vo that incorporates a major portion of 
the second spherical harmonic and leads to a closed solution with no secular variations of order k. The perturbations of 
this non-Keplerian intermediary orbit are derived here by the von-Zeipel modification of the method of Delaunay. The 
secular variations are carried up to orders k? and k’, and the periodic variations to orders k and k’/k. A computational 
‘summary is included; the effect of the third spherical harmonic is appended. 


1. Introduction. Until the advent of the arti- 
ficial satellites, the restricted Particle-Spheroid 
problem has held the distinction of being the 
‘simplest of the unsolved problems of celestial 
mechanics. Under the assumption of axial and 
equatorial symmetry, the gravitational potential 
of an oblate planet may be approximated by 


| V=-—1/r+ 2kP,(sin 6)/r? + Rk’ P,(sin 6)/r®>. (1) 


Here ¢ is the radius vector, 6 the declination, P» 
and P, the Legendre polynomials; the planetary 
constants k and k’ are assumed to be small 
quantities of the first and the second orders 
respectively. 

Since there is no exact solution of the problem 
of motion in spherical coordinates, attempts have 
been made to find a V» that leads to an exact 
solution and that approximates V better than 
the term —1I/r. The author found such an 
intermediary orbit (Garfinkel 1958). In its sim- 
plicity it approaches the ellipse; in its accuracy 
it far surpasses the latter. A definitive solution 
of the problem must include the periodic varia- 
tions of O(k) and the secular variations up to 
O(k?). The author first solved the problem by the 
method of Poisson, based on the Taylor series 
expansion, in powers of k, of the equations of the 
variations of the constants. In this paper the 


results are re-derived by the more elegant and 
powerful method of Delaunay. The latter is 
based on a canonical transformation that re- 
moves the periodic terms from the hamiltonian. 
The solution includes as a special case the recent 
results of Brouwer, obtained by the application 
of the same method to the elliptic approximation 
Vo => — I/r. 

2. The intermediary orbit. It has been shown in 
the earlier paper that 


YVo= —1i/r 
++ 6k ¢1(sin? 6 — ¢y)/2r? + ¢3/r + ca/r? | (2) 


is the unique potential with the two properties: 
(1) It achieves the Hamilton-Jacobi separability 
in spherical coordinates and leads to a closed 
solution in terms of elliptic functions. (2) It 
preserves the gross features of V; i.e. the equa- 
torial and the axial symmetry, the singularity at 
the origin, the vanishing at infinity, and the 
order relation V + 1/r = O(R). 

T. E. Sterne (1958) has treated the case 


¢3 = 0; the author has chosen the case c, = 0, 
Vo = —u/r + 3kci (sin? 0 — co) /7?, 3) 
bw = 1 — 6ke3, 3 


which minimizes the number of elliptic functions 
in the solution. With this choice of Vo as a first 
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approximation to V, the disposable parameters 
C1, C2, and cz; have been determined so as to re- 
move the secular variations of O(k). The prin- 
cipal results of that paper will be briefly sum- 
marized here. 

The equations of the orbit defined by (3) are: 


E-esnkE=nt+oa, 


r =a(I — ecos#), 
v I+e ] E 
tan - = tan =; 
2 I—e 2 
u = (vu + w), (4) 


sin @ = sinzsn u, 


¢ — 2 = (Az cos t/dj) 
x f du/(1 — sin? isn? x). 
0 


Here a, e, 1, ¢, w, 2 are constants analogous to the 
usual elliptic elements, and 


Pp rar a(i oe e”), ve = By 
€ = 6ke./up, 
AP = 1+ e(c2 + cos? 2), (5) 


ho? = I + e(C2 — sin? 72), 


the modulus « of the Jacobian elliptic function sn 
being given by 


Ko (E/N SIM eer. (6) 


(In the original paper, with the results carried 
only to O(k), the factor \;? in (6) was omitted.) 
In terms of the circular functions the last three 
of the equations of the orbit become: 


Y= (1+ gn)\@+~), = 
v= P+ desin ay t---, 
sin 6 = sinzsin y, (7) 


¢ = 2 + tan“ (cos 7 tan y) Ad gga) cine se, 


where the constants go; and g3. are related below 
to the fundamental frequencies of the motion. 

The three fundamental angular frequencies 
have been given by 


ny = i, 
Ny, = mari /2K, (8) 


= 
| 


K 
= (nd2 cos 1/K) i du/(I — sin?z sn? x), 
0 
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where K(x) is the complete elliptical integral of 
the first kind, with the series expansion in x’, 


2 4 . 
=2(1+£4% 4...) @ 
The quantities 


Zor = Mo/m — I, 


< 
10, 
232 = N3/n2 — I, ( / 


in view of (8), (9), (5), and (6) can be approm 
mated to O(R*) by 


ee Caer ce 64 


ll 


4e(3 Cos? 4 = 1 +263) 
+ drel—5 + 8c. — Bez! 

+ (18 — 24c2) cos?4 — 21 cos*z] 
—}(e/dx”) cos 4 


— ;:2(5 cosz7 — 3 cos? 2). 


(Ir 


£32 


The orbit can be characterized as a pseudo 
ellipse ; the motion in 7 is identical with that in < 
Keplerian ellipse with Vo = —y/r. The vari 
ables E, v, W are the analogues of the eccentric 
anomaly, the true anomaly, and the argument 0 
latitude. If c; = 0, then x = 0, go: = g32 = 0 
1 = Ae = uw = I, and the orbit degenerates int 
an ellipse. 

If the c;’s assume the values 


I ; ; | 
a= ? C2 = cos? 1, (12 
c; = (3.cos?i — 1) V1 — &/4p?, 
the secular variations of O(k) are zero. Witl 
the abbreviations 
6k 
SS P =4(1 — 3cos?2), | 
CN ap tos ) (13 


Q = <sin’1, 


the quantities \1, A», M, Ze1, Zs2, Which involve th 
c;’s explicitly, become on the substitution fron 


(12): 

AZ = I + 2€ cos? 1, Ke = 
w= 1+ 6RPVI — &/p°, 

to O(k), and | 


I + €cos 22, 


ga = Fi (5 cos? z — 1) 


| 
| 
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2 
ar (— 5 + 26 cos?z — 53 cos‘z), (15) 


2 


€ € : 
| SS —— — , —S4{( — , 3 4 
£32 5 cost + Fe ( 5 cost +19 cos’ 2), 


to O(k?). The non-Keplerian orbit defined by the 
above choice of the c;’s will be designated as the 
“intermediary” and will be taken as the un- 
perturbed orbit in the Delaunay theory. The 
en problem, defined by k’ = 0, will be solved 
first; the perturbations arising from the fourth 
spherical harmonic will be calculated in section 
10. 

| 3. The old Hamiltonian. The Delaunay vari- 
ables for the intermediary orbit are: 


b= f nit +, 


| L = Vua, 


16 
( GaINi=#, emo 
k=, 


| JEL NS CSO ie 


Except for the appearance of the factor \2, which 
is nearly unity, they are identical in form with 


those for a Keplerian ellipse with Vo = —p/r. 
|The Hamiltonian is 
F = 2/21? + R, (17) 


where —y?/2L? is the undisturbed total energy, 
and the disturbing function R is to be expressed 
in terms of L, G, H; 1, g, h. From (r) and (3) 


R = —2k(P,» (sin 0)/r° 
| — 3¢, (sin? @ — c2)/2r? — 3¢3/r]. 


(18) 


In the analysis the c,;’s must be shown explicitly 
until all the necessary differentiations with 
‘respect to a, e, 7 have been performed. With the 
|Bpbreviations (13), the replacement, of sin @ by 
sin z sin y, followed by the replacement of y by 
‘b+ 4eQsin 2¥ +---, in accord with (7), con- 
verts the last expression into 


m= Rit R.+:::, 
re cos 2y — P) oe 


i 
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where the subscripts of R indicate the order of 
magnitude, and y has been written for y. A 
standard form for R as a double Fourier series in 
v and y can be derived by the replacement of 
b/r by 1 + ecosv. 

In the application of the Delaunay method it 
will be necessary to construct the secular and the 
long-period parts of R and of its derivatives Rz, 
Rg. Let — be the slowly-varying quantity 


f£=y—v2 e010 + g, (20) 


that is the analogue of the argument of the 
pericenter go in elliptic motion. If now a function 
f(, &) has the Fourier expansion 


f= ¥ Agcos (i + 232), (21) 


then its secular part f and its long-period part f* 
are defined by 


f =A, f* = X Aojcos 2jé, 7 #0. (22) 
0 


The Fourier coefficients in (22) are given by the 
integrals 


I 20 2a 

ai f fdldé, 

it 27 27 

af i f cos 2jédldé, 7 #0. 


The change of variables /,  — v, € may be carried 
out with the aid of the law of areas, 


= (nr?/G)dv, 


A Ne See 
(23) 


(24) 


which is rigorously valid only in the unperturbed 
orbit. However, the resulting error is of a higher 
order, so that, to O(kf), (23) becomes 


it Qr Qr 
caf, J) oeseyanas, 


a. hy 1) (nr? f/G) cos 2jédudé, 
27 J0 0 
j- A 0. 


Avo = 


(25) 


Aoj — 


It follows immediately that the principal parts 
of Aoo and A; are respectively the constant term 
and the coefficient of cos 27£, 7 # 0, in the ex- 
pansion of ur? f/G into a Fourier series in v and &. 
If such an expansion is available, the integration 
(25) is unnecessary, for A oo and Ao; can be found 
by direct inspection. The following results will be 
useful : 


(= Gree 
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Ri=F ~~ 4P + op(Q — o/2) 


r r 
(26) + cp'/Vr—@], (7 
cos 2y = 0. pe 2 
R, = op (1 — &)!Qap(cip — 1). 
The secular part of R now follows by the substitu- 
tion from (26) into (19): In terms of the Delaunay variables, 
ain ac LG eee a6" ( Ele ) 66;G°L 
ag “| - z ¥ 2 Gr? a 2u ate Go? ar we : Ca 
als oR’c H \2 k 
2 TGs (G7 = "4) (: =e 4 OR), i 
with R* = 0; (31) 
2 = 1+ 6Ral(co —1)G2 + HG], (29) i.e., the disturbing function has no long-period 
part. 


in view of (5), (13), and (16). With f = Rin (25) 
observe that in the expression 7*R the terms 
involving y and r are of the form cos 2¥/r and 
cos 4¥/r. Consequently, in the expansion 


PR = Yi Axe, Cos (Riv + 2kow) 
= > Ai; cos (iv + 27&) 


k. = O, 1, 2, and if ky # othen |k:| = 0, 1. Since 
4 = ky + 2ke and j = ko, the inequality 7 ¥ 0 
implies 7 ¥ 0 and Ao; = 0, so that 


(30) 


The derivatives of Ri with respect to L and G 
which will be needed later, are obtained with the 
aid of the relations 


Ri = Rid, + Ret Re = Rea + Rta, Ga 


the literal subscripts indicating the arguments ot 
partial differentiation. The quantities Ru, R., R 
to O(k), and the Poisson brackets @z, er, @a, 1¢, 
given in the earlier paper, pp. 91-92, yield on 
the substitution into (32): 


PR G, eG, (P — Q 2y) = ee ) ae es 
(PRit/ hee 30 cos 2y) | e Ae cos v — e cos 2u mS! 


2 2 
— Qsin 2y (4 ap =) sin v + 3e sin 2v + = sin 30| 
+ $c:p(4Q sin? y — c2)(— 3e + 2 cosy + e cos 2v) 
+ c:pQ sin 2¥ (4 sin v + esin 2v) + csp? (cos De 22" (33-1) 


P 


(P-Rie/G)i = lea — 30 cos ay) E se (1 + #e?) cosv + e cos 2v + - cos 30| 


2 


+ Qsin 24| (4 ain c)sine + 3esin 2v + Ssin 30| — ecos*isin? y(1 + ecosv) 
— $c:p(4Q0 sin? y — c2)(e + 2cosv + ecos 2v) — cipQ sin 2¥ (4 sin + esin 22) 


Since the terms cos (2¥ — 2v) cancel out in 
(33.1) and (33.2), 


(Rit*)1 = (Rie*)1 = 0; (34) 


i.e. the derivatives of the disturbing function 
with respect to L and G have no long-period 
parts of O(k). It must be noted, however, that 


+ cipe cos? 7 sin? Y — csp” cos of (33-2) 


the integration of a long-period term with respect 
to the time is accompanied by a reduction of the 
order of magnitude by unity. For this reason it 
is necessary to carry the analysis to O(k?) in 
order to obtain the long-period term to O(k) in 
the final solution. Then the approximation y, = I 
of the earlier paper must be replaced by 
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> = 1 + ga + g20 cos 2y. The term go: yields 
| Jong-period parts of O(R’) : 


(7? Ri1*/G)»2 — €£210(G/L) cos 2é, 
(P-Rig*/G)»2 = + €£210 cos Pe, 


jwith € = y —v and go = O(e). The term goi0 
‘cos 2& contributes merely short-period terms of 
| O(R*), which do not concern us here. 

| The treatment of go: as a constant is justified 
| by the circumstance that the dependence of go: 
/on G contributes only a pure periodic term of 
O(R?) to angle y. Indeed, this contribution is 


(35) 


= (gai)e@ f dGdy, 


inasmuch as, to the second order, 


SC Ag A yd(8G). 


Since 6G is purely periodic, it follows that the 
mixed secular term in the element g is cancelled 
in the perturbation of the position coordinates. 

4. A method of Delaunay. In the old Hamil- 
‘tonian, F(L, G, H;1, g), the variable kh is absent 
in consequence of the assumed axial symmetry; 
the time ¢ is absent, having been incorporated 
into /. The conservation of the axial component 
| of the angular momentum, H = const., and the 
conservation of energy, F = const., follow im- 
mediately from Hamilton’s canonical equations. 
The problem is thus reduced to two degrees of 
‘freedom. We seek a transformation to new vari- 
ables L’, G’, I’, g’ such that /’, g’ are absent in the 

new Hamiltonian F’. Then L’ and G/ will also be 
constants of the motion, and /’, g’ linear functions 
of the time. The secular variations can then be 
‘derived from F’, and the periodic variations from 
‘the transformation function S. 

In his Lunar Theory, Delaunay used a suc- 
‘cession of canonical transformations to remove, 
‘one by one, the periodic terms of the Hamilton- 

ian. The feasibility of a single transformation to 
accomplish the same purpose was suggested by 
von Zeipel (1916). Its first application to arti- 
ficial satellites was made by Brouwer, who used 
the elliptic orbit, c; = 0, as a first approximation. 
For any choice of the c,’s, the generating function 
S(L’, G’;1, g) of the desired canonical trans- 
formation satisfies the relations 
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L=S;, G=S,, 
Uk — Sine Ca = Sq. 


Since R is of O(k), S must differ from the identity 
transformation, 


So= Ll + Ge + HL, (37) 


also by a quantity of O(k), and it must include 
terms S; of short period, as well as terms S,* of 
long period. Accordingly, let 


Sa CG ee-p at Bebo) ores, 
Si = OCR’), 
and, up to O(k?*), 
‘Sys 1S ES eet a erie 
Slee ig 51, ay aa 
Since S* contains / and g only in the combination 
£ = g20(1) + 2g, 
Sir = Lo10151," = OR), 
Sor*® = go1d1S29* = O(R®), 


(36) 


(38) 


(39) 


(40) 


thus justifying the omission of S2:* in (39). The 
new Hamiltonian must satisfy the relation 


MCE Gt, 2) + 3S: = FL G): (41) 
In view of (36) and S; = 0, (41) becomes 
Sih Oye bre) Gs (42) 
Now let 
P= Fo Fy + Fos +5 Fe = OCR’), Ga 
Peelig anh Oop lou tata ty 
where we identify 
Fo fel, Fy = Ri, Fy ={Rx 44) 


The substitution from (43) and (39) into (42), 
followed by an expansion of the left-hand 
member into a Taylor series about the ‘‘mixed”’ 
variables L’, G’, 1, g, yields, up to O(R?): 


Fe Fa al: 
= Fy(L’) oF For(Siu ce Poe hi S21) 
+3 Ports + Pit PirSu 


deePialSuetyoig Jateehe, (45) 
In view of (40) and the relations 
For = —”, Fort = 3n/L, (46) 


separation of (45) by the orders of magnitude 
yields 
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PF = fig pefar??; 

Fy = Fy, — 1Su, 

Fol = Fo + FyrSu + FieSiy (47) 


+ (3n/2L’) Si? + (Fie — g210)Sig* 


= nSo1, 


with the arguments L, G replaced by L’, G’. The 
secular part of (47.2) is 


Fi = Fy 


since $;,; is purely periodic. 
determined from (47.2) as 


= R, (48) 


Now 5S; can be 


I — 
Sy = =a (Ry = Ryydl 


n 


= if (7?Ri/G)dv — Ril/n. (49) 
The secular part of (47.3) is 
Fi! = Fy + = R, + 4, (50) 


where @ is the secular part of ©, defined by 
® = FiSu + FigSiy + Gn/2L Si. (80) 


The remaining terms in (47.3) contribute nothing 
to ®, since 6 and /;¢ have no long-period parts, 
while S*,, is purely long-periodic. Therefore the 
long-period part of (47.3) is 

) 0 = (Fig — g210)Siy* + #*, (52) 
where p* is the long-period part of &. The ex- 
pression in parentheses equals —mgzo, since 
0 = n, and since the vanishing of the first-order 


secular variations on the intermediary implies 
Fig = o. Now let 


i cf pas aj sin 27&, 


b* = 5° A; cos 27€, (53) 


3ky 


Pee 8 se H 
1 TAC Ra ee 


= $ne COs 1(cip—1), 


with the abbreviations 


M, = cp(Q <= C2/2), M;= csp?/VI ae Ce, 
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= Trea 144 0- S) (aa, 


= neVi —e(—P+3M.,+2Ms3), 


sae SVD 36,G” 
i’ = S| -344(5- a )( 


= [ne E cos*41—1+4¢\p (p- a} 


JBf y+ poe 
Gd» 


HH? I 
yaa 19 — 
a) + f 61G’2(1 202) 
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where £ = gow + g with L, G replaced by L’, G 
wherever they occur as arguments of v. Then ‘Se 
can be determined from ' 
a; = A;/2jnga1, (54) 
A; being the coefficient of cos 2j(y — v) in 
Fourier expansion of nr’@/G in y and ». 
The new Hamiltonian F’ now furnishes 
the secular variations through the Hamilton 
equations 


/ = rh PN ia = By h! = — Fy’, 
while the periodic variations follow from the 


equations 


ie eons Li + S11, 

G= 2 = G" ate Sar 42 Sige , 

Hie aa wo i Fe 
l = 1+ Sin + Si 
g’ deeb =g+ Sig + Sig, 

W = Sw =h + Sin + Sie, 


involving the derivatives of the generating 
function S. Explicit expressions for these varia- 
tions will be developed in the following sections. 
In the process the values of the c;’s will be re- 
determined. ia 

5. The redetermination of the c;’s. The dis- 
posable parameters are to be determined so as to 
remove thé secular variations of O(k). The latter 
are obtained by the differentiation of Fi’, which 
is given by R; in (28.1) with L, G replaced by 
L’, G’. The derivatives of \2, which appear it 
this calculation, being of O(k) and multiplied by 
k, may contribute only terms of O(k?). Accord- 
ingly, with dp» fixed, 


cae 
(1 —2¢2) +4¢3 —j—- | 
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and with a, e, z redefined so that 


L'=WVya, G’ = L'vi — 2, 
| (59) 
HT = )2G' cos 1. 


For the vanishing of the first-order variations in 
(57) it is necessary and sufficient that the c,’s 
assume the values in (12) with a, e, 7 redefined 
as in (59). For c; = 0 (57) reduces to the elliptic 


values: 

i’ = 4ne(3 cos? i — 1)V1 — e, 
_ G1’ = 3ne(5 cos?z — 1), (60) 
i hy’ = —4necosi. 


| 6. The new Hamiltonian. The new Hamilton- 
ian is 
r 


a F’ = 2/2L°+Rit+R,+%, (61) 


‘with R, and R, furnished by (28) and L, G re- 
placed by L’, G’. To evaluate & from the defini- 
tion (51), expressions for (2iz)1, (Rie)1, Sir, Sty 
i required. The first two are furnished by (ar 
the last two will be now derived from S; in (49). 
‘Mirst, with the aid of (19.1) and (27.1): 


PRG = &| @ cos 2 — ap)(2 ) 


+ cip (20 sin? y — 2) + cat (2) (62) 


R,/n = = «GL — 4P + ap(Q = C2/2) 
=e C3p?/VI +i e |. 


| 

| 

hen the replacement by I + e cos v of p/r 
in KGS. qi is followed by a conversion to a Fourier 
series in yw and v. Prior to the integration of the 
resulting expression it is necessary to replace L 
and G by L’, G’ wherever they occur as argu- 
ents, implicitly and explicitly. In order to avoid 
“mixed” variables here, as well as in other 
algebraic operations throughout the analysis, a 
total replacement of all the old variables by the 
new will be made whenever Z and G must be 
replaced by L’, G’. This practice is justified by 
the fact that the resulting error is of an order 


higher than our tolerance. The integration of . 


(62.1) is further facilitated by the use of the 
approximation y = v+ g, with g constant to 
O(k), and of the relation 
f 5M —édv=E=1+esinE, (63) 
0 ‘ 
which is rigorously valid only in the unperturbed 
motion. Again, the resulting errors are of a 


"ae 
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higher order, so that, to O(k), 


S,:=«G | —3P(v—I+e sin v) 
+40 sin 2vy+esin (2¥—v)+ sin ev-+s)| 
—4¢,pQ sin 2¥+cip (o- 2) (v—l) 


+ (csp2/V1 —e) esin x}. (64) 


Now the differentiation of S;in (49) with respect 
to 1, followed by the use of (19.1) and (62.2), 
yields S\,; the differentiation of S; in (64) with 
respect to g, which enters only through y = v + g, 
yields Si,: 


Su == (Ri — Bi) 
3 
oe #14030 cos 2y—P) (2 ) + 4P(1—e?*)8? 
1 “2 b \’ 
+ 3c1p(4Q sin? ¥ — ¢2) ( ”) 
(65) 
marci (O =" 69/2) (r=. 67)! 
+ Csp* (2 mai ec )| 
Sig = «GQ E cos (2W — ») 
+ (1 — ¢1p) cos 2¥ + = C08 (ay + | 
The last term of (51) leads to 
Cromer at is ae Sar 
a oe aan Rt Rv 
= 5 Ri- Ra, (66) 


which is easily evaluated from R; in (19) with 
the aid of (26), and from R, in (27). Therefore & 
is the constant term in the Fourier expansion in 
v and yw of the expression 

LN Rise RigSiy ose R Re. (67) 
G (RinSu + RiaSig) onl tr). 7 
The multiplication of the double Fourier series 
for Riz, Rie, Siz, Sig is perfectly straightforward ; 
after some laborious algebra, the second order 
portion F,’ of the new Hamiltonian is found as 
the sum F,! = R, + &, with LZ and G replaced 
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by 25 'G: 

py) ee ats {— + x + 5x0? 
YA 32 LG" 5 4 


+9? (10 — 24% — 18x”) +y*(35 + 36x + 5x”) 
+6¢7°G"4[3 — 8¢.+8c?+ y?(8c2—2) —y* ] 


+240G”[—1+2¢.+y?(2—6c¢,)—y*}}, (68) 
with the abbreviations 
x= G/L’, y = H/G@d2. (69) 
We have omitted in (68) the portion 
18k? 
(My + M;)(3M2 + M; — 2P)x, (70) 


L3G" 


inasmuch as it contributes nothing to the per- 
turbations of the intermediary. Indeed, from 
(58), (12), and (13) 


M, = — M;= P, 


the first-order derivatives of the entire expression 
must also vanish. 

7. The secular variations. Having removed the 
secular variations of O(k), we proceed to calculate 
those of O(k?). Note that the quantity po, enter- 
ing Fy’ through R, in (28), contributes nothing 
to the intermediary orbit since the derivatives 
of A» are to be multiplied by the factor cyp — 1 
which vanishes by (12). Furthermore, the deriva- 
tives of \» are of O(k); consequently A» entering 
F,’ through FR», contributes only terms of O(k*). 
Accordingly, it is sufficient for our purpose to 
differentiate F,’ in (68), treating \» as a constant, 
and then substitute the values of the c,’s 
from (12): 


. i ~ 
a oe foe Bor Wei Se Bye > 


+ 3°(138 — 96% — 90x") 
+ y4(— 129 + 144” + 25x”) |, 


EY 


I 
go = >>> n€lL— 101 + 24x 4+ 25x? 
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’ I ‘i 
ho’ = po ae 23 + 12x + ox? 
+ 9° (43 = 30% eae 


with x = VI — é?,y = cos2, « = 6k/G* 7 = ae 

For c; = 0 (68) leads to the elliptic values, 

(i ms =e en + 6x + 2 ; 
2)0= ona 15 + 16x + 25x , 


+ 37(30 = 96% — 907) 


¢ 
aS y4 (105 + 144x a 25x") |, { 
se hi 
(g2’)o = [— 35 + 24% + 25x 
a7 (73) 
+ y? (go O20 a 126x?) 
(iul)o = — Hols — 12x - 924 


+ 9°(35 ++ 36x 4 See 


in agreement with Brouwer. 
8. Variations of Short period. From (50) the 
variations arising from Sj are ; 


6L 
él 


Si, 6G = Sib 
Sunny, 02 =) Sas, 


sf = 0, ‘ 
bh = — Sin. (74) 


The derivatives Si, and S,, are found directly by 
the substitution of the values of the c,’s from (12) 
into the expressions for Si; and Sy, in (65); 
Siz, and Sig are obtained by the partial differ- 
entiation of S; in (49). In view of the relations 
n = wL-, dl = (nr*/G)dv, and the vanishing of 
R,, Re on the intermediary, the latter operation 
yields: 


; 
7 


1 


Sin = BS yl la it (??-R1/G)dv, 


wn 


See f (PRa/G)db. (75) 


Tey (72) The values of the c,’s can now be substituted 
5 ; from (12) into the expressions for-S1, Riz, Rig in 
+ y'(462 — 192% — 126x?) (64), (33), (35). With the retention of the long- 
+ 9 ( (407% 200% 4 A570 uF, period terms of O(k?): 
Si = eG{3P[2v — 3E + 1 — esinv] + 4Qe[sin (2y — v) + sin (2¥ + 2) I, | 
2 
r-R1r/G = — (eG/eL) {P| «- I +-'ViI — e — 42) cos wi -+ 7008 3v — 2e (Zv —-@— :)| 


0 


2 


2 
| - 4008 (2y — 3v) + 2eg21 cos (2W — 2v) + (— 1 + Ye?) cos (2p — v) 
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+ (— 1 + 48e") cos (2W + v) — 4ecos (2y + 2v) — = cos (ay + |, 


rRe/G = — (é/e) |p t Sie 2 — #e?) cosu — S608 30 (76) 
+ BI cos (2b — 3u) — 2ega1 cos (2¥ — 2) + (1 + Je”) cos (2p — v) + 4e cos 2p 


Ba, a e*) cos (2p + 2) + 4e cos (2p + 2u) + _ cos (2¥ + 30) | 


2 


ap S cos? i [eos 9 — 4cos (2W — v) — 4cos (2y + o)}]. 


| The integration in (75) is facilitated by the use of the relation (63) and the approximations y = v + g 

for the short-period terms and wy — v = gow + g for the long-period terms. The resulting error in 6/ 

and 6g arising from all sources will be of O(k?). Finally, Sin is obtained by the differentiation of Sj in 
| (64) with respect to H, which enters only through z in P(z) and Q(z), with 77 = G~' cot 7, followed 
by the replacement of c:p by unity. With the abbreviation 


|| 


| ri 3 La nee) eerie) er dima (eore)(e — x)4(1 > x) (77) 


| the results are collected below: 


{ 2 2 
Blog L = @/(r - |p Je+ (: - r= 4 - ©) cose — S cos 30] 


| 


+ QO cos 2y E + (1 + 2e?) cosv + e cos 27 + 3008 »|}, 
«eQ [cos (2 — v) + 3cos (29 + »v) ], 


| 6 log G 
él = (eVI — e/e) |p| (- 1+ v1 — & — 2”) sin» +— 5 sin 30 + & sin F] 


Si ae C3) p ecin (2y — 20) -- (— 1 + $e”) sin (2y — 2) 


i E 


+ (— 4+ xe?) sin (2y + v) — esin (2p + 22) — =sin (ay + so) ||, (78) 


| —————s e ; 

| 6g = (€/e) P| a — Vi — .— 3e%) sine — 7p sin 30| 

| 

| + | - sin (2y — 3v) — esin (2y — 2v) + (1 + fe?) sin (2¥ — v) + 2e sin 2y 


+ (4 + +e’) sin (29 + v) + esin (2¥+2v) + sin (2y + 30) | 


2 


+ - cos? 7 [sinv — # sin (2y — v) — Fsin (2¥ + 2) J, 
dh = dee cosi[— sinv + 4#sin (2y — v) + ¢sin (2y +») ]. 


In addition to non-singular long-period terms 9. The long-period term. From (56) the varia- 
in (78), there also appear in the solution terms tions of long-period arising from S\* are, to 
with the singular factor (5 cos? i — 1)~, arising O(R), 

} from S*;. Such terms will be considered in the 65L*=0, dG*¥=S,,*,  5H*=0, 


next section. = Sir) 08 = — Sia, Oh* = — Sin (79) 
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To obtain S:* from (53) and (54) it is necessary 
first to evaluate the long-period part of ® from 
the definition (51). This evaluation is equivalent 
to the determination of the coefficients A; of 
cos 27(W — v) in the products of double Fourier 
series. In the result, 


@* = A, cos 2& 
I . 
ire néG'e? sin? 1 (I-15 cos? 2) cos 2& 
A,= As= +++ =0, (80) 


curiously, all the terms involving the c¢,’s ex- 


Syt= eG’ (r= 2) (197) (Ea) Aa tana 


k G? gre oe eee 
= F6GAN Tat Ge ae CaN ia 


Finally, the differentiation of Si*, followed by 


the restoration of ¢, e, y, £in the right-hand mem- . 


bers, and the use of the abbreviations 


Fa (i a ae 
sg, = rere 
Big Tae te oy ‘ 
Spa (83) 
yields: 
6 log G* = ee’a cos 2&, 
6l* = — e(1 — e)8 asin 26, 
dg* = € i(: = “) O71 «| sin 2&, 
6h* = — ee’y sin 2¢. (84) 


It is noteworthy that the coefficients of the long- 
period terms, being completely independent of 
the c,’s, are identical with what they would be 
in the elliptic approximation c; = 0; indeed, 
(82) and (84) are in perfect agreement with 
Brouwer. Incidentally, (84.2) corrects the 
_ author’s previously published erroneous value 
of the long-period term in the mean anomaly 
[1958 ]. 

The anomalistic and the draconic frequencies, 
n, and m2, become equal at 7* = tan-!2 = 63°4, 
giving rise to the phenomenon of “‘resonance,”’ 
with singularities appearing as poles of orders 


the secular part of AR appears as 


/ 


r= e9)8(1 + 82) Ps (coss) = 
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plicitly and all the higher harmonics of 2 have 
vanished. The divisor mgo1, occurring in (54), is 
seen from (15) to be 


(81) 


to O(k). This expression is identical with the 
secular rate of the advance of the argument of 
the pericenter go’ in a perturbed ellipse, c; = 0. 
That the trigonometric argument = gow + gis 
the ‘analogue of go follows from the equality of 
their secular rates and, to O(1), from the equality 


of their initial values. Now Si* can be written | | 


ngo1 = 4en(5 cos?i — 1) = Go’ 


:) | sin 2(ga0-+g). (83) 


one and two. The long-period is 


P = r/ngon = 4n/ne(5 cos?z — 1), (85) 


with 


min P = 1/6k (86) 


corresponding to @ = 1, e = 0, 7 = O. For the 
earth with k = 0.545 X 10~-* and the canonical 
unit of time 806.8 sec., min P = 8.9 days. 

10. The effect of the fourth harmonic. The con- 
tribution of the fourth spherical harmonic in the 
potential to the disturbing function is 


AR = — k'P, (sin 6)/r°®. (87) 
From the definition 
P, (sin 6) = 3(3 — 30sin26 + 35 sin‘) (88) 


and the relation sin @ = sinz sin w there follows: 
P, (sin 6) = 2Pz (cos 2) 


+ 3, sin? 7 (6 — 7 sin? z) cos 2p ; 


+ oe sin‘ 7 cos 4y. (89) 
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Upon the replacement of P4 Ce 0) and (p/r)5 Pa 
their secular parts, 


P, (sin 0) = 3P, (cos 1), 
rod SE 2 
Gin = a —2(1+%), (0) 
‘ 
Fee GC Hw Ht 
ee Fal cai 35)(s a 3° Ge an 6S =), (91) 


Ee ee 


| 
| 
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_and the contribution to the new Hamiltonian is 
| AF = AR, (92) 


with L, G replaced by L’, G’. From AF” are de- 
‘rived the secular variations, which can be written 
in the form: 


Al! = — dne'x(1 — x) (3 — 309? + 359%), 
Ag’ = gen<'[— 21 + 2707? — 38544 
+ x2(9 — 126y? + 189y*) ], 
Bal = 4ne'(5 — 3x*)y(7y? — 3), (93) 


‘with the use of the abbreviations x = V1 — e?, 
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y = cos7, and 


e = 15k'/8p*. (94) 


In view of (87) and (89), the discussion of 

section 3 makes it clear that the long-period part 
of AR arises only from the expression cos 2/7’. 
Indeed, 
AR* = ine’Ge’ sin? z (1 — 7 cos? 2) cos 2, (95) 
with & = y —v = gow + g. This expression for 
AR* appears as an additional term in the right- 
hand member of (52); the contribution to S,* is 
therefore 


/ 
LAS,* = =(£) G’e? sin? 7 (1 — 7 cos” 2) (5 cos?z — 1) sin 2 


is 4 


| 4 I GP? 

: = 4 €1p Gi 1 aoarg L2 I 
with the constants e: and p defined by 

6. = G4= 6k, p= ¢/e = 5k'/96k?. (97) 


| From (96) are derived the long-period variations. 
| With the use of the abbreviations 


ea ey h(t 74") (Sy? —~ 1), 


\ 


[8 = dae = 97 + 89" — 97, 8) 
| my = B'/4, 

| they can be written, analogously to (84), 

lB log G* = epe’a’ cos 2, 

| 6l* = — ep(1 — e”)3?q’ sin 2, 

| 2 

| ég* = ep I(: + ) ea ee | sin 2£, (99) 
| 6h* = — epe?y’ sin 2é. 


|The short-period variations arising from the 
fourth harmonic are of O(k?) and can be neglected 
| in our analysis. 

| 11. Computational procedure. The secular 
variations of O(k?) can be incorporated into the 
orbit equations (7) by replacing 1, go1, g32 by n’, 
Zo’, Zs2' defined as follows: 


n=ntl, + al, 
= ga + Go’ + Ag’)/n, 
B32 = Bae + (hg! + Ah’)/n. 


The substitution from (15), (72), (93) into (100) 
yields the values that appear in the summary 


below. 


(100) 


ll 


"a 


dak Jake 


We eee 
~s +725 \(sa-*) sin 2 (gv +g), (96) 


It is convenient to apply the periodic varia- 
tions in L,G,/,¢ directly to the coordinates, 
with the aid of the first-order differential 


relations: 
éy = r,6a + r.6e + 1/01, 
oy = 6v + dg = v,6e + v,6l + dg 
= (5g + d1/V1 — e) + v,6e (101) 
— [(r — oVi — &)/ejedl/Vi — & 
with the partial derivatives 
hen Sk, Fe SS GLO, 
ry = aesin v/Vi — &, 
ve = (2 + ecosy) sinv/(I — e”), 
v, = avi — &/P. (102) 


lt is now clear that the singularity in 6/ and ég 
in (78) at e = 0 has been removed. The varia- 
tions 6a, 6e, 62 are related to 6Z and 6G by 


6 log a = 26 log L, 
be = — [(1 — e)/e]d log (G/L), (103) 
6t = cot76 log G, 


which follows from the definitions (16). Finally, 
explicit expressions are derived for 67, 6y, 62 by 
the substitution from (78), (84), (99) into (ror) 
and (103), followed by the substitution from 
(103) and (102) into (101). The results appear in 
the summary below, with the terms arising from 
Si* marked by an asterisk. 
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SUMMARY 
A. The data (for a satellite of the earth) 
1) the geophysical constants: 
R= 6378.4 km., mo = 1-23943 X 10 rad-/sec: 
Jo = 2k = 1.083 X 107, 
i 


2) the elements a, e, 1, a, w, Q 


— 1.2 KX 107° 


B. Auxiliary constants: 


=i (654, 


Fe 
= Ao 
I1+<*x 


8 
| 
4 
| 
x 
< 
| 


a, = $(1 <= x) ay 
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C. The secular effects: 
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P= 7 — oy7)" 0 = 7a 
b=1+ J4/J? 
bb) = 1 — 4b 
m = qybo — xb(59 — 1)" 
a = 4mQ 
y = yL—debo + 3065 — 1) ] 
B= V9. 

GAP = 50 1d) 
p= ax’, B= 1+ 3x ee 
€ = 3Jo/up, m= wPa*?no. 


+ y4(— 129 + 144% +:25x") — 72p(1 — x7) (3. = 3097 ease ee 


+ yi(— 815 + 360x + 45x") + 240[— 21 + 27077 — 385y* + 47(9 — 1267) en sea 


n’=njt+ a [— 33 + 16% + 25x? + y?(138 — 96x — 90x") 
gor = Ze(5y? — 1) + rr {— 131 + 24x + 25x? + y?(618 — 192% — 126x”) 
7 é 
Bea = Ores ale 53 + tax + owt 97(157 136% — 5x?) + 240(5 — 3x") 7 ee 


D. The equations of the orbit 
l=n't Se Ws 
Eee hessinme 


v Ite iy 
tan - = (aN =a 
2 f= 2 


= a(I — ecos £), 


= (1 gor) +e), 


‘ 


Bee |- Petes sm eneosen 
+ | £ cos (2p —"0) + 4 c0s 2 
+ © cos (24 + |} 
OL a ca| Seos (2y — 2v) + ecos (2p — v) 


e2 
+ - cos zu, 
rv =rt+ or + br*, 
6Q = Zee cosi[— sinv + 4sin (2p — »v) 
+ § sin (2y + »)], 


6Q* = — ee’y sin (2y — 29), 


OF = QO + 60 + 6n*, 


by = e{P (2a: sin v + ae sin 20) 
+ QLe sin (2¥ — v) + = sin 2p }} 
— cosz 6Q, 
oy* = « E sin (2y — 2v) — 2easin (2y — v) 


ea . i 
— — sin 2 
2 ’ 


Z 


Y=et oy + sy*, 
Care 
ein 21 [cos (2y — v) 

+ 3 cos (2 + v)], 

62* = deme? sin 27 cos (2 — 20), 

v=itéet i, 
sin 6’ = sinz’ siny’, 

¢’ = Q + tan (cos 7’ tan p’) + gao’p’. 

The position of the satellite is furnished by the 
perturbed spherical coordinates 7’, 6’, ¢’. The 
equatorial radius R is the unit of distance; for 
computational convenience the time ¢ is here 
expressed in seconds; the term $x’ sin 2y, be- 
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‘longing to the unperturbed orbit in (7), has been 
‘incorporated into 6y. The calculation of tangents 
could be avoided by the use of the alternate 
| relations 


avi — é 


‘ a 
sin EZ, cosv = ; (cos E — e) 


¢ = Y + cos (cos y’/cos 6’) + 230'p’. 


| It is noteworthy that the solution appears in 

| closed form; there is a singularity at 7 = 7,, but 
_ mone at ¢é = Oori = 0. 

; The constants J;, introduced by King-Hele 

_and Merson (1959), are the coefficients of the 

| spherical harmonics in the potential V. They are 
related to the Jeffreys constants and to our k and 

| k’ by 

i” 

| 


J, = — $;D = Tels 


12. The validity of the solution. The solution is 
valid under the following restrictions: 1) e < I, 
i 2) |t| <7, where 7 is a bound beyond which the 
omission of the secular variations of O(k*) can 
no longer be justified, 3) |i — i,| > w, where 
the bound w is the half-width of the “‘resonance”’ 
| region in the neighborhood of the singularity 
at. 2,. In particular, w will be so defined that 
| outside the region 
[y*| < ec = ef’, (104) 
| and consequently all the long-period terms can 
_ be treated like quantities of the first order. If 
|¢-1i,| <1, the formulae of the preceding 
| section furnish the principal part of the bound on 


ey], 


|6y*| < e|2ea + e6]., (105) 
With the use of the abbreviations 
= Aue (106) 


A= 4|5cos?s—1| > |i —i,|, 


the dominant terms of the asymptotic series 
for a and 8 as cos? iz — # are 


= 7a, |p| = C/4d’, 


In order to assure the inequality (104) ihe its 
therefore sufficient that 


(107) 


|a| 


Ase SG VC + 4Cp*)/p. (108) 
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If it is now assumed that 
b = O(1), (109) 


then (108) can be approximated by A/e > VC/2p, 
so that w can be chosen as 
fee VC (110) 
i= Pe 
In particular, for an artificial satellite of the 
earth the relation (109) is certainly satisfied ; 
furthermore, the launching conditions generally 
imply p=1-+e. Then (110) becomes ap- 
proximately 


w = e|b|?/20(1 + e) (111) 


with max w = 1.°4|b|? corresponding to e = 1. 
The singularity at 2, dissappears when w = 0. 
It is noteworthy that this occurs when e = 0 
and also when 6 = 0; the latter is characteristic 
of the Vinti potential (1959), with J, + Jo? = o. 

In a separate paper it will be shown how a non- 
singular solution valid for all 2 can be constructed 
by the retention of g’ in the new Hamiltonian. 

13. Comparisons and checks. Checks of the 
theory are provided by a comparison with the 
Brouwer solution of the main problem. Three 
such checks have been carried out, and perfect 
agreement obtained as far as the analysis went: 

1) The amplitudes of the long-period varia- 
tions are identical with those of an ellipse, c; = 0. 

2) Secular variations of O(k) and of O(k?) 
reduce to those of an ellipse when the three 
disposable parameters c; are set equal to zero. 

3) Since the perturbed intermediary must be 
identical with the perturbed ellipse, the three 
natural frequencies must be the same in both 
orbits. These frequencies are: 


ni =Vant+i/ +l, 
ni = : = (1 + gor) (m1 + gi’ + G2’), 
ie = b = (i — £32) No" a he +- Jigt : (112) 


i ae from the facts that the bound of 0 
and the total energy 7’ must be the same in the 
two orbits, it can be shown that the constants 
L’, G’ in the two orbit are correlated by 


AG’/Gy Same lee 
L'/Lo = pli — ¥A(ePG/L) + LAF,’ ], (113) 
where AG’ = G’ — Gy’, etc.; and the subscript 


o refers to Brouwer’s perturbed ellipse. Since 
cost = H/)\2G’, and n = p?/L*, it follows from 
the definitions of \», P, «, and from (68) that 
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AG er) ah : AR = — k3P3 (in 0)/r4 
ake « 4 ais P3 (sin 6) = 3(5 sin? @ — 3 sin 8) 
Ae = — 4&P, 42 sin? i — $ sin 7) sin 
A(ePG/L) =3;eV1—e (1-4 cos? 1-+7 cos‘), AB — 3sin*zsin 3p 
=O 
~ = — (ePG/L)o If ks = O(k?), only the long-period term need be 
0 


+ 2 exte—69 + 1599 


= «1 + OF). (114) 


With the aid of these results and the expressions 


for go1, Zs, U, g’, h’, it has been established that 
An; = Any = An3 = O. Indeed, the anomalistic, 
the draconic, and the sidereal frequencies of 
the two orbits are respectively equal, as they 
must be. 

The relative size of the variations in the two 


orbits is of interest. Not only do the secular ° 


variations of O(k) vanish on the intermediary, 
but the short-period variations 6Z and 6G vanish 
with the eccentricity e. The analysis of section 4 
shows that the expressions for F,’, F3', etc., 
involve the powers of 6L and 6G. It is therefore 
clear that the smallness of the latter quantities 
generally assures the smallness of secular varia- 
tions of O(k?) and higher. Indeed, for a typical 
orbit with e = 0, 7 = 45° the solution of the 
main problem yields the values 


= fetal ae LT is 
bie FET (G 


eee Fo 
which are only 3% and 6% respectively of the 
corresponding values 


17V2 


(he')o = — 192 


en 


on the ellipse c; = 0. The relative smallness of 
secular variations of O(k*) on our intermediary 
‘should guarantee the superior accuracy in the 
calculated position. This advantage is achieved 
without any sacrifice in the ease of orbit cal- 
culation. 


14. The effect of the third harmonic. Our as- 
sumption of equatorial symmetry can be removed 
by the inclusion of the odd spherical harmonics 
in the potential V. The effect of the third 
harmonic will be calculated here by the method 
of section 10: 


considered, which arises only from the expression 


sin y/r*: 


AR* = +3nk3;G~e sin 7 (5 cos?7 — 1) sin &, 
AS\* = + 3(k3/k2)G’e sin 7 cos &, ‘ 
= 3(k3/ke)(1/G? — 1/L”) 
Xx (i. — 1/6?) cose 
kop = 2k, E=pP-—vSgvt+¢ 


With the use of the abbreviation 
Ge k3/2R2G? 


the resulting variations in the Delaunay elements 
can be written: 


6 log G* = + ese sinzsin &, 


e6l* = + e3(1 — e)? sin7z cos &, 
sin 4 
6g == a( S24 —e 


bh* = 


— ee cot 7 cos &. 


With the end of (101), (102), and (103) of 
section II, there follows: 


ér* = + e3G? siny, 


éy* = + al — + (2cosv — esin?v) sin] 
sin 7 
X cos &, 
67* = + ee cos7 sin &, 


C= 4 4 Ott: 


The remaining singularities in 62* and éy* a 
4 = 0 are removed in the variations of the 
coordinates 6 and ¢: 


sin 6’ = sin7z’ (sin y + cos péy*) 
= sinz’ sinw + e3cos €cosy 
X [e + (2 cosv — esin? v) sin? 7] 
é A 
cos? 66¢* = — 3 E (sin? v — sin? y) — cos | 


X sin 27 cos ~ 
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formulae, which led to the discovery of the mixed 


secular term. 


‘The variations 67, 67, 6 sin 6, dy can be readily 
incorporated into the computational scheme of 
‘section II. 
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|Dr. Dirk Brouwer, the Director of Yale Uni- 
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jencouragement received throughout this work. 
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|the conscientious checking of the mathematical 
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| THE MOTION OF A CLOSE EARTH SATELLITE 
: 


Abstract. In the present paper perturbations of six orbital elements of a close earth satellite moving in the gravita- 
| tional field of the earth without air-resistance are derived as functions of mean orbital elements and time. No assumptions 
) are made about the order of magnitude of eccentricity and inclination. However, it is assumed that the density distribu- 
| tion of the earth is symmetrical with respect to the axis of rotation, that the coefficient of the second harmonic of the 
| potential is a small quantity of the first order and that those of the third and the fourth harmonics are of the second order. 
| The results include periodic perturbations of the first order and secular perturbations up to the second order. 

| However, the solutions have some singularities for an orbit whose eccentricity or inclination is smaller than a quantity 
| of the first order, and this case is treated in a different way. 


in the expression of the semi-major axis. 

1. The disturbing function. In the present paper 
it is assumed that air-drag is absent and that 
‘the gravitational field of the earth is axially sym- 
‘metric. Under these assumptions the gravita- 
‘tional potential of the earth at a point where 
‘the geocentric distance and the latitude are 7 
and 6, respectively, is expanded into the series 
of spherical harmonics, 


p GMI) Ast int) 
= , {r+ 2(4 ~sints 


A 
+ 22(Ssin's — 3) sins 
r\2 2 
4+ 24(43 + <sin® 9 
— sin 
SO 


— © sin’ 20) te of, (1) 


where G is the gravitational constant and M is 
the mass of the earth. The second and fourth 
terms in the expression of the potential are due 
to the oblateness of the earth, and the third term 
is due to the asymmetry with respect to the equa- 


By using Delaunay’ s canonical elements a theorem is proved that there are no long-periodic terms of the first order 


torial plane. A» is taken to be of the first order 
of small quantities, and A; and Az, are to be of 
the second order. The coefficients of higher har- 
monics may be of the third order of small quan- 
tities or less. (O’Keefe, Eckels and Squires 1959, 
Kozai 1959). 

The purpose of the present author is to derive 
the periodic perturbations of the first order and 
secular perturbations up to the second order. 
Therefore terms of higher order than the fifth in 
the potential series may be neglected. 

As the satellite is always on an ellipse, whose 
position and shape are variable, it is convenient 
to express 7 and 6 in (1) by elliptical elements of 
the satellite by the following relations: 


a(t — é?) 
1+ ecosv’ 


. (2) 


sin 6 = sinzsin (v+ a), 

where @ is the semi-major axis, e is the eccen- 

tricity, 7 is the inclination to the equator, w is 

the argument of perigee and v is the true anomaly. 
The disturbing function due to the oblateness 

of the earth is then 
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| 
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A Sa I ieee ey 
(2) {+ — Esints +5 sin? cos 2(0 + «)} 
a\r Awe 2 2 


4 
) (3 SI ee 3) sity Fe) = ® sin? Z sin 3.845 »)| sin ¢ 


I 
—~sin?z + S sin' i + sin? z (2 = 5 sin’ i) cos 2(v + w) 


ee > sint 7 cos 4(0 2 a | . @) 


The true anomaly, 9, is easily transformed to the mean anomaly, M, which is a linear function of time 
in non-perturbed motion, by the differential equation 


dv a 
Wit (4) 


r/a and v appearing in the disturbing function of (3) are then functions of e and M only, and are 
periodic with respect to M. Therefore, R is also a periodic function of M and w. In R, terms depending 
neither on M nor on w are called secular, terms depending on w but not on M are long-periodic, 
and terms depending on JM are short-periodic. 

As the long-periodic perturbations originate from terms of the second order in R, we must retain 
secular terms and long-periodic terms up to the second order. However, for short-periodic terms 
we need only terms of the first order. 

Using the following relations (Tisserand 1889) : 


(Gas  (G) se 
Ta 27 Jo Ne 
a 3 
(*) cos 2v = 0, 
ie 
PENI. Aakers 
(¢) cosv = e(1 — e?)—5/2, 
aNe Gane 
{ (7) cos 37 = (*) sin 37 = 0, (5 


anan(r4 de), 


aaa 
218 
eee 
w 
a 
=) 
iN) 
a 
I 


2/18 

~~ <z 
f = D 
x18] =a 
—< a 
ll I 


5: 
(*) cos 20 = oe — e712, 
r 4 


GND a \? GP. 
= sit 20e— COS AU — nen Sine — nO) 
: B 4 = 4 , 


we can pick up only necessary terms by the above criterion as follows: 
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R, = uS(4 — = sin? i) (1 — &)-3/2, 


uo 
Ry = ne (1 +3e), 


A 
$2 sin (Sin Bien :) NC ie Sin. 18 (6) 


eee ee LO 
+ otsints (3 — Ssine i) er — @)- "cos 20} 


Azfa\ I I : 
ia om (2) I(: — 5 sin’ a Vf _ (2) (I = ys] + = sin® i cos 2(u + »|, 


| where Rj, Ro, R; and Ry are first-order secular, second-order secular, long-periodic, and short-periodic 
parts of the disturbing function, respectively. 

2. Perturbations of short period. The differential equations representing variations of orbital 
elements are: 


| eas 
| dt nadM’ 
ik det 6 ok VI —e&dR 
| dt nae dM nae dw’ 
dw _ COs 1 ok Ni- @orR 
dt naV1 — sini % nave de’ 
| (7) 
j di cos 2 aR 
dt na2Vi — e sini Oo ’ 


dQ I OR 


} Tae neni — e sini Oi’ 
| OM | b= OR +2 OR 
| dt a nae oe na 0a’ 


where x is related to a by na = GM. 

| To derive short-periodic perturbations of the first order, after replacing R by Ry in (7), one may 
| regard a, n, e, 7 and w on the right-hand sides of equations (7) to be constant. However, n, appearing 
as the first term of the last equation without any factor, is variable, but it is a known function of 
‘time, after obtaining the expression of the semi-major axis. True anomaly, v, may be regarded also 
‘to be a known function of time on the right-hand side, and one can transform the independent variable 
from t to v by 


dt Ne I 
d= dM = — =(2) dv. 8 
dM 5 Pape (8) 
Then, for example, the short-periodic perturbations of the inclination are obtained by, 
; cos 7 eNO 4 
diy = n’a?(1 — e?) sin a(t (<) dw ae (9) 
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where the integrand can be expressed by a finite trigonometric series and is, of course, integrable 
analytically within the necessary accuracy. The results for the six elements are as follows: 


CA 


9 3 . 
da, = 2 (1 _ 8 sine) {(£) — (I - a} ++ (¢) sin? 2 cos 20 +0) |, 
Ge NB 2 r r 


JS pe 9 3 3 
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BUNS 2 r 2\ 7 
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Senee yl « e) si ee ee 
seit i 5 sin’ t ; I 4 e)sinv + > sin 2v = sin 3v 


eee i: Jie eine 
— sin? z + 5 epee e?+ sin (v + 20) + TE sin zsin (v — 2w) 


St ee ae |cos 210 +6) + ecos (v + 2) + = e.cos Ge + 2a), 


Sees een C2. oe ane 
+ gsin®7 sin (4v + 2w) + 7g sin’ 7 sin (5u + 2w) |, 


Ndi. 7 sin 22 {cos 2(v + w) + ecos (v + 2w) + = cos (3v + 2a) ; 


( _ S sina) s sin 2(v-+w) +— ae sin? 7 — = (: _ = sin? i) e| sin (3v-+ 2w) (10) 


| 


dQ, = fo — M + esinw,— = sin 2(v+w) — < sin (v + 2w) — esin (39 + 22) ; 
dM, = t= -(:-3 ne 1-©) ‘ide cian eae 
F sin? z F si 5 sin 7p Sin 3 
+sints|2(1 + Se) sin @ 4 24) — Sein @ — 24) - Z 2 (30 + 20) 
sin ; i sin (v ++ 2a) a6 sin) (W920) 5 “Smee sin (3v + 2w 
A 
3 P 
—3e sin (4v + 2w) — — g sin (5u + 22) |, 
- where 
= a(I — e). 
As the mean values of cos jv (j = I, 2, ---) with respect to M do not vanish, but 
5 ae 7 - 2 
cos jv = (——} (i + jV1 — e), ce 


mean values of these short-periodic perturbations are not zero, except for those of a. 


values with respect to MV are: 


Their mean 


—— 2, 
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aa rm sin? 2 COS 2U COS 2, 
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A, >; I I —@ I Je ee a 
dw, = Pp sin 3 la pope COS 20 Nt = COs 2 COS: 20). Sin 2u, 


6é 6 
a= aly Ao 
a, = — 5 Pe sin 22 COS 2¥ COS 2, (12) 
I A» 
MQ, = a4 P COS 2 COS 27 SiN 2w, 
e 
A ie Pons 
2 Ls £ 5 
dM, = — a ioe é sin? | 5 + ea 608 20 Sin 2w. 


Therefore, de, — de,, dw, — dors, di, — dis, dM, — dM, and dQ, — dQ, must be the short-periodic 


perturbations whose mean values with respect to the mean anomaly are zero. 


Expressions of the mean anomaly and the argument of perigee are rather complicated, so it is 
better to combine the four elements a, e, w and MM into the radius vector 7 and the argument of 
latitude L = v + o by the following relations: 


dr é : r da 
ee pa — cos v de, 
2 
ao = Fr = eau + sin o(1 +5) 2 ae 
Putting 
As 
da = — aoe (1 — 3 sin’ i) vi — @& + da,, ea—nee 
A» 
dw = dw, — a ae sin? 7 sin 2w, 
3 Ae APIS 
dM = dM, ee VI — e’ sin? 7 sin 20, 


_ the deviations of the radius vector and argument of latitude of the satellite, from those computed 
_ by mean orbital elements are obtained as follows: 


Be Bas 3 sine‘) eae ane ii soa uc 
re (x — Sines I at I ISSO 6 a cers gee ~cos 2(v + w), 
A 


ee [(2 a S sin’ i) ee Aen v) 
} (13) 


V2 ( en: (r= Wr =) sin 20} 


Ss a Dale Pen aky 1 Ogee 24 od | 
6 sin i) ésin (v + 2w) € 7p sin i) sin 2(v + w) 6 cos # sin (3v + 2w) |]. 
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The secular perturbations of the first order are easily derived by putting R = R; in (7), as 


ae 
at SO aae 2 = Sina ab 
= Aa, _ 
Q = Q — pt COS 2 
M = M, + iit, 


= 
| 


= mo + Go no (1 — Sin Wire 


where wo, Qo and My are mean values at the epoch,-that is, the ‘aitial values, from which periodic 
perturbations have been subtracted. is the unperturbed mean motion, which is related to the 
unperturbed semi-major axis @o by nc’ao’ = GM. 

It is more convenient to adopt as a mean value of the semi-major axis not do but 


A 
d= 09 {1 - 2 (1 — Ssint) r—él, 
p 2 


so that the following relation holds: 


Vor GM x = A ~ 3sinei) vi =ah. (14) 
And to derive the expressions of (13) this value of @ has been already adopted as a mean semi- 
major axis. 

3. A theorem on non-existence of long-periodic terms in the semi-major axis. As is well known, Poisson’ 
proved a theorem of non-existence of secular perturbations of the semi-major axis of a planetary 
orbit. However, there are some differences between theories.of planetary and satellite motions. 
Because the perihelion and node of an ordinary planet moves around the sun with the period of 
some ten thousand years, one may expand e¢ sin w, e cos w, 7 sin Q and 7 cos ® into power series of 
time in the planetary theory, so that terms depending only on w and Q must be regarded to be secular, 
not long-periodic. From this point of view Poisson proved the theorem. Now we observe the motion’ 
of an earth satellite for a long interval of time, in which the line of apside makes several revolutions. 
So the corresponding theorem in the satellite motion is that there are no long-periodic terms in the” 
expression of the semi-major axis. 

Let us transform variables from: Keple r’s elements to Delaunay’s canonical ones: 


Ib = Vua, G = Vua(i — e), lal = Vua (1 — e*) cos 1, (as) 
i= M, [Gaby [pe == . 
where up = GM. 
These variables must satisfy the following canonical equations: 
dt dF dG Vans) eae uae 
dt ol’ db 1de'" dt OR? | 
(16) | 
ws WP \ tae OR ah ek : 
GEV anh Sau en mee OG?) wads oH’ . 
where 
2 
Pa aR 
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Since F does not depend on h, there is an integral H = constant. Variations of the semi-major 
axis are represented by: 


Cok eORs 
Eno ae tore 


where y?/2L?, Ri, R2 and R; are omitted in F because they do not depend on /. 
Long-periodic perturbations of the first order of Z will come from long-periodic terms of the 
second order of 0R4/d/, which may be expanded into Taylor’s series, 


ay (eee ores aR PRs ee 
al ( al ) + (Sa )az + ( an ) a+ (See) ac + (Fs decane (17) 


where terms of higher than the third order are neglected and in parentheses /, L, g and Gare replaced 
by their respective mean values, / = at + 1, L = Lo, g = gt + go and G = G; dl, dL, dG and dg 
are deviations of instantaneous values from their respective mean values and are of the first order. 


nw and g are secular motions of the mean longitude and the longitude of perigee, and are supposed 


to be known. 

In equation (17) one must consider only long-periodic terms of the second order, and after inte- 
grating with respect to time only such terms of the first order, omiting other terms. The first term 
of (17), 0R./dl, cannot have any long-periodic terms of any order, so must be omitted. As R, depends 


_on ¢ only through M within the accuracy of the first order, the following relations hold: 


O7R, a I ORs if O?R, ¢ I of) 
fe). | alee , 
S(Ge)a-3(), f oe) at = 2 Re) 

0loG ni \ 0G dldg n\ dg 

Integrating 0R,/d/ in (17) with respect to time by parts using the relations (18), one has 
OR, 
i ig 
i ORs a) ORg OR, 
al (ar) + (Se) + (Ze) a0 + (31) | 
I OR, dL OR, dl is OR, dG OR, dg ; ) 
Ailes dt + al Meet) tae) a ie Woe aid } a oe 


After integration the first part contains only terms of the second order, so can be omitted. The 
second part is transformed to 


w= E6008) (28) +28) fo) 
: + Ge) Ge) + Ge) Ge) 
tf an( SS) at os 


dL 


ll 


ll 


: | 
( 
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By using the relation nL’ = const., this equation is written as 


OR, 
dL = L $f (| ) dLdt 


as S [RaL] 2 al (Ry) Sat 
= 3 crt —% f cry (%) a 


- strana ef (Fa | 


). (21) 


ahiae Biba 
“VRE [Rees 27 


Since the terms are all of the second order after integration, one can conclude that there are no long- 
periodic perturbations of the first order in the expression of the semi-major axis. Of course there 
are no secular terms in the semi-major axis. 

4. Secular perturbations of the second order and long-periodic perturbations. Let E; be one of the six 
orbital elements of a satellite, and express its variation by the differential equation 


oe 


dE; 

\ 

A function f; may be expanded into a power series of deviations from the mean orbital elements as 
oe a | 

=) +(e) am + -. (23), 


Secular perturbations of the second order and long-periodic perturbations of the first order will come 

from the terms of the second order on the right-hand side of (23). However, the complete first order 
expression of dE; is not yet known except for da. By integrating equation (23) by parts as in the pre~ 
vious section, there holds, ; 


z= fi (fildt + pe [FijdE;] — X mf Fa as re ye oa | 
where : 
y= {®) dt 


_ If £; is the inclination, 7, for example, iat 


f cos 4 oR : 
a = . . ie” f 
n@VI — esin 7 Iw | 


In this case 0f;/0E; has neither long- periodic nor secular terms as far as the terms of the first order 
are concerned. As it is easily proved that 7 and e cannot contain any secular terms, >); /ijdE; can 
be omitted from the equation of the inclination, because this is of the second order and has no” 
secular terms. 

Therefore all terms on the right- hand side of (23) for the inclination are known and long- periodic 
terms can be picked up by using the following relations (Tisserand 1889) : 
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3 
(*) cosjv = 0, (j > 1). 
Then the expression of long-periodic perturbations of the inclination are derived as: 


ecosisinw. (25) 


yee Alias e* Sin 27 14-15 sin?z A, 18-21 sin?z 3 As 
diy; = a; = 


ce = chess 
p® 8(4-5 sin? 2) 6 LE Oe, paras vase 
' Now there is an integral, Va(1 — e”) cosz = const., and it has been already proved that there are 
no long-periodic perturbations of the first order in the semi-major axis, so de:, the long-periodic 
perturbation of the eccentricity, is 


ie" sit 


ns Pe cost m 
ee oe 1 e.sin? 7 14-15 sin?z A, 18-21 sin?z i De neti 
= de, + Pees cin?) 5 ri Ae 7 cos 2w + a sinzsinw. (26) 


However, the same principle cannot be applied to obtain the expressions of the node and the 
argument of perigee, because in these cases Of,/07, Of;/de and df;/da have secular terms. But 
fortunately complete expressions of diz, de and da have already been derived. And since it is also 
proved that Fj,.dw and F;ydM do not include secular terms, the following expressions are derived: 


a Az. Z Ax 3 e 2 i i(8 S42 4 ‘| 
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where 7 and é are mean values of the inclination and the eccentricity over all the periods, and 
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It is very difficult to derive the long-periodic perturbations of the first order of the mean ano 
because the long-periodic perturbations of the second order of the semi-major axis are not kn 
However, usually the atmosphere of the earth changes the mean motion of the actual satelli 
much that any long-periodic variations of the mean anomaly cannot be detected from observa 
with good accuracy. 

The results are then: 
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> € and 7 are mean values with respect to 
», and wo, Qo and Mo are initial values from 
1 periodic perturbations have been sub- 
od. On the right-hand sides of these expres- 
one must always use constant mean values 
e, 1,» — ot, Q — Qt and M — it. 

The case when the inclination or eccentricity 
y small. In the previous sections it is as- 
d that 4-5 sin? z, 7 and e are not very small. 
s very small, periodic perturbations of M 
» become very large. However, in the case 
> short-periodic perturbations the difficulty 
ranish, if elements are transformed from 
»torandv+o. 

r the long-periodic perturbations the vari- 
must be transformed from e and w to 


—& = e cosa, 
n = —esinw. 


2 variables must satisfy 
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> » has the same value as in the previous 
on. 
e solutions of these equations are 
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d wo are constants of integration. 

ey is much smaller than 343;/4aA.2, which is 
2 first order, the argument of perigee cannot 
around the earth completely, but oscillates 
id the value 90° as follows: 
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The variation of the eccentricity is expressed by 


uesly ee aie 
é=>—~sinz+ eosin w. 


- 4 aA2 


In this case dQ; = di; = 0, because they have a 
factor e. 

When the inclination is very small, instead of 
w, the longitude of perigee, r = w + Qis adopted. 
The longitude of perigee moves secularly as 


aT = wt oF T 0, 
where 


é A» 
J > Nn. 
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The variations of 7 sin 2 and 7 cos 2 are derived 
as in the previous case as 
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where 
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If 7) is much smaller than 3¢43/4A op, then 
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Abstract. Sections 1-6 give the solution of the main problem for a spheroidal earth with the potential limited to the 
principal term and the second harmonic which contains the small factor kz. The solution is developed in powers of ke in) 
canonical variables by a method which is basically the same as that used in treating a different problem by von Zeipel/ 
(1916). The periodic terms are divided in two classes: the short-period terms contain the mean anomaly in their argu-_ 
ments; the arguments of the long-period terms are multiples of the mean argument of the perigee. 


| 
The periodic terms, both of short and long period, are developed to O(k2) ; the secular motions are obtained to O(k»”). | 
The results are obtained in closed form; no series developments in eccentricity or inclination arise. The solution does not 
apply to orbits near the critical inclination, 63°4, but is otherwise valid for any eccentricity < 1 and any inclination. 
Section 7 gives the long-period terms and the additions to the secular motions caused by the fourth harmonic in the ~ 
potential; section 8 gives the contributions by the third and fifth harmonics; section 9 contains formulas for computation. 


1. The equations of motion. The equations of 
motion of a small mass attracted by a spheroid 
are 


Py OU Wy OU Veen ee 
ad? dx’ df@ ey’, dP az’ 
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waft yee es ea 
= as A (i = 2:sin*)p) 


+ (1 — ro sin? ¢ + © sin* ) +--+ 


The equatorial plane of the spheroid is taken as 
the xy plane; 6 is the latitude, and if M is the 
mass of the spheroid and k the Gaussian con- 
stant, uw = RM. 

The main problem of artificial satellite theory 
may be considered the problem with k, = o. Let 
I be the instantaneous inclination of the orbital 
plane with the equatorial plane, -g the’ distance 
of the pericenter from the ascending node, and f 
the true anomaly. Then 


sin 6 = sin sin (g + 15 
2 sin? 8 = sin? J [1 — cos (2g + 2f)], 
and the disturbing function can be written 
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Let a, e be the osculating semi-major axis and 
eccentricity, respectively. The Delaunay vari- 


ables are then 


L= (ua)}, /=mean anomaly, 
G=L(1-—e?)}, g=argument of the pericenter, 
H=G-cos I, h=longitude of ascending node. 


With these variables the equations become 
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The parts of F not exhibited in terms of the 
Delaunay variables may be expanded in Fourier 
series as follows: 

a? T3 ie) 
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The coefficients P;, Q; are power series in the 
eccentricity ¢; the lowest power of e is e? in P;;_ 
el¥2I in Q;, 7 +0. Since e? = 1 — G?/L*, the 
derivatives of a function ¥(e) with respect to L™ 


or G may be obtained by ; 
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The functions y that will arise are functions of 
a/r and f. Well-known formulas of | elliptic 
motion are 


2(2) =F c0s f (4) 
LI? 
of (242) sins: (5) 


Some important properties of the Hamiltonian 
F may be noted: 


(a) the independent variable ¢ is not explicitly 
present in F, hence the integral F = constant 
exists ; 

(b) the variable # is not present in F; 

(c) the coefficient Qo in a2 is zero. 


The value of this coefficient is the constant term 
of the expansion in terms of the mean anomaly 
of a*r~* cos 2f: hence 


I T Qs 
Qo -:{ ps 008 2 fal. 


By making use of the integral of areas in the form 
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dl = Giadf 
this may be written 
ele a 
Q=—G : 7 cos 2fdf 
pee (.~ 
=12f (1 + ecos f) cos 2fdf 
10; 


This derivation may be generalized to show that 
the constant term of (a?/r?) cosqf (p and g 
positive integers), expanded in terms of the 
mean anomaly, is zero if g > p — 22 0. 


2. Outline of the method of solution. Consider a 
transformation from the variables L, G, H, 1, g, h 
_ to new variables L’, G’, H’, I’, g’, h’, with the aid 
of a determining function S(L’, G’, H’, l, g, h). 
Then, if 
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as as os 
L= al’ a Og’ H= on 
es 0s ? as i! OS 
Mora MOCK Omen Te: 
the equations in the new variables will be 
dL' oF* dl One 
Chemo sowed atid COL 
CGO ea oF* 
dt re og’? dt cs OG (7) 
CHEE GIR all! rele 
Gr mons di am No lr es 


with 
FEC OG AH 31" 9! =) 
= (NG Hele) 1(8) 


The dashes in the places for h’ and h are used to 
indicate the absence of these variables. 

The problem would be completely solved if a 
determining function were found such that F* 
is a function of the variables L’, G’, H’ only. The 
differential equations (7) show that then L’, 
G’, H’ are constants, while /’, g’, h’ are linear 
functions of the time. Substitution of this solu- 
tion of the primed quantities into (6) then yields 
expressions from which the original variables 
may be obtained in terms of L’, G’, H’, I’, g’, h’ 
and therefore in terms of ¢ and the constants of 
integration. 

In the present problem it is more convenient to 
choose the determining function S in such a 
manner that /’ is not present in F*, while g’ is 
permitted to appear. If this is accomplished, L’ 
and H’ will be constants, and the system is 
essentially reduced to one of one degree of 
freedom: 


dG’ art dg’ _ 
Cie moe dL 


oFf* 
dG’ 


in which L’ and H’ are present as constants. 
After this system is solved, /’, h’ are obtained by 
quadratures from 


dl! OF* 
ai Ob). dio 


Preferably a canonical transformation by the 
choice of a suitable determining function leads 
to an equivalent result. 
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3. First-order solution. The determining function 
is obtained by a method used by H. v. Zeipel 
(1916) in a qualitative study of the motions of 
minor planets. In this procedure a development 
in powers of the small parameter » is introduced. 
If F, is written for R&, the original Hamiltonian 
is 
P= Fo bly 

in which the subscript denotes the power of k» 
present as a factor. Mo is a function of ZL only. 
Let also 


iS == Scrap Dalsts Oo oaey 
FY = Fy* - Fy a eh 
It will be convenient to choose 
So=L14+ Gg t+ Hh. 


If then the expressions for L, G, H,l’, g’, h’ 
given by (6) are substituted into 


FL, G, Hel, 8) =F" (4) Ga 
there results 
n (#2) «(A 
Ntpalge @ Gf = - ; —) 
ey fe GU < ; -). 


Expanding everywhere by Taylor’s theorem to 
the second power of ko, the result is 


aha eSs 
OF o OS» I 0? Fo OS? ears 
TOL! al +350 ( 5) 
=P HEM EDE. Ge. JEM L, g, =) 
OF AS) OF: AS 
Aso, aae noe 
a Fy* Se NE AGRE Gi. is ~, g, =H) 

OF 6S: 

dg aq + Pe (CLs Cw cha ae 


Parts of corresponding order in ky on both 
sides yield: 
order 0, 
Fo(L’) 


leh Nc (9) 
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order I, 

oF OS 

aL Sern ke 
order 2, 
dF) OS» ay a a(S 2 en OF, OS, 
dL’ al 2. 02 No! OL’ dl 

OF, 0S; OF y* 0S; 

aG! ag == Ht ag IG” (11) ’ 


The expansion could be carried on indefinitely, 


but for current practical requirements considera-' 


tion of parts beyond the second order in ke ap- 
pears to be unnecessary. Put 


Py = Fis ae Isis 


in which F\, is the part independent of /, Fi, the 
part dependent on 7. In the terminology of 


(10) 


oe 


planetary theory /;, would be the secular part | : 


of the disturbing function, /}, the periodic part. 
Hence, if for the sake of brevity the notation 


Fike ae Sak aga: HT” . 
= es PR eae B= 3 ae (12) 
is used, 
uke 
Fis = L3G3 ’ { i 
za te 
fy = (Aoi + Bor). 


Now equation (10) is split up into two equations 


wk 2 


ar = L®’ (Aoi te Bo»), , 
(13) 
4k 
Fi* = nate A, 
LG" 
in which use is made of 
OFo _ ge 
yy ope IGS 


By integration S; may be obtained in the form 
of an infinite series as 


= 


wh 


= Fs 


Eee ae - P; sin jl 


+ BE Qsin (e+ in|. a 


No constants of integration are required, since 
only partial derivatives of S; with respect to 
L', G’, H’, 1, g, h will be needed. 
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| vantageous to obtain a closed expression for S. 
| The method used at the end of Section 2 for 
| obtaining the integral of functions of the type JL ie 

| ar-” cos gf can be used to obtain a oG [; sin (2g + 2f) 


foa-f(S-G)a 
oidl = 2 @ + > sin (2g + f) + @sin ce +3/)| 
L 
= Li? + esin f], 


For the following developments it will be ad- | Rte f s Hee onal 


tn 


Hence 


2h } 
Sy= a= 1+esinf)+B | $sin Cg 27) =p < sin Ch OD eae <sin (2g + pI, (15) 
| 
ios ah 
| ‘ie 7B | cos eg 4 27) +e cos (2¢ + f) + cos ce +3/)| (16) 
Introduce 
H Uke 
Voy = L'*? 
then 
acer 
5 = Jb, eee a 
: Pa NN O87 8 2 Res - 
ee la a) ig )emerten|, an 
We 2S 


ey f 4 v Bal 2 - 3) [cos (2g + 2f) + ecos (2g + f) + = cos cet+3p]| (18) 
Aside (19) 


S, was obtained as a function of G’, H, e, f, g. Hence L’ is present through e and f only. To obtain 
| 0S;/0L’ it is convenient to obtain fret 

Lo: 2k a 

° = 4 | « + ecos f) ot sin ‘| + B cos (2g + 2f)(1 + ecos f) a 


+ sin Age aia sin (2g + s/f 
in which use was made of 
cos (2g + 2f) + 5 cos (2g + f) + 5 cos (2g + 3f) = cos (2g + 2f)(1 + ecos f). 
With 
ra) 2 CG 
(+ ecos 5! = (SE + 2) sin f 
there results 


OS; whe aG”? Ff 
Boe = Gi \4(S5+¢ + :) sin f 


9 2 a? G' : 
ae |(-She-S+1)snce+n+(Sht$ ae £) sin @e +3|I. 


| 


anal 
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With the aid of the relations (3) there follows at once 


jee ae SZ \(SS a )si C 37) 
: ral v3 Ge Pa Aan ut 4 4G? 


2 i 
(-S35-f41)sin cet n+ (GR+o+3 +) sin ee +3), (20) 


rie 
Lt 
Bes: ih 
£ > 8 aG’ oa 
ey ye ee ee Eten es 
=f BE £4 3H) (EF +241)sing + (3-25) 
a? G”? a é WiGle ee le Nae 
x |(-S&- $+ 1)sin e+ +(Goa+$+5) sin Ce + 30] 
Ibe 2 aot Serr 15 
+nen{(-3+4 a) tes - 57) 
x E sin (2g + 2f) + sin Gey a qsin (2g + s/f, (25) 
Bee: 
rec 0H . 


=~ 3nZa[p— 14 esin s —Fsin (e+ 2p) —Ssin @e +) — sin ee +30 |Z. G2) 
In the calculation of the coordinates, / and g are not needed separately but only the sum: 
il? L” as) Lt K I a Er a? G” 4 sa 
b+g=7 +2 4 > “(e =) {( 1437 )(S5 42 +1)sin + (3 - 4G3 
a? G” i aG? i : 
‘i te Pi aa r) sin meg ee +244) sin ce +3|| 
L" Bo ER EE referee 2: 
* 1 Za {(— 2 +, 2 E) (f= re eine) Via (2- 3) 


x E sin (2g + 2f) + £sin ee +P + £sin (og + s/f. (eas 


Since (I — e?)-! — (1 — e)-+ is divisible by e?, the first part of the difference (J + g) — ie + g’) has 
eas a factor and not as a divisor. 4 


4. The second-order terms. For use in the evaluation of F2* it will be convenient to put 


po = 7 cos (2g 27) a €:cos (22 Fiat )eep 5008 (2g + 3/)| (24) 


so that 
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| Also, define 71, 72 by 


Bees 3 a it 
cee eee 2 I~ 3 Ge as) 
25 
ies i |G page SG . GPIL 
eee be AES ay 2 & | cos Cae “| § mu 2S | cos (2g + 3/), 
| which are easily obtained with the aid of the Peas ey and (5). Then 
oF 4h 
| a on ale 6(Aci — Boz) pains L? Laie + Bn) (26) 
OF 4k Jet JEP 
a= ~ Gs | Fant Br) +35 Fe — os) | @7) 


If the terms dependent on / are not needed to the second order in ks, the only part of interest in (11) 


will be 
I 0?Fo oy OF yy OS, a7 S| A 
[SS/ a) ) ob’ al ° 0G dg eee 
_ where the subscript s designates the parts independent of /. The contribution to the right-hand 
member that has 0F,*/dg as a factor vanishes because /\* is a function of L’, G’, H’ only. 
The following products are easily obtained: 


1 (95:)P _ 3 ath 
20L"\ al 2 BM 


(Ao; + Bo)? 
OF yp 0S p®Re 
Ai eh) aa oe 


| OF aS: pike bee 


— 6(Ao, + Boo)? + Fa (Aor + Bos) (Ani + Bra)| 


L” 


(Gl 
— o1p2 + o2p2) — = (ABrips + B*rop2) 
L ; 


AGL dg L"0 7 Gn B 

| OF 29S: p®ke? i 

| eee ‘Elana :) | 
| and hence, 

pk? 9 G” 

| Fy* = a {4° (- MENS aa nn) 


G” G’ Lt 
+ AB | — Qoi02 + Fr (c1T2 ae 7271) es jee Ta oc 3 Ga P2 
G’ 
[P(-28+& sours — & raps) 


JEANIE ee 
Dray al a (— 301p2 + 30292) — 3 ios |h . (28) 


The parts independent of / may be found by eet integrals of the same type as those which 
occurred in the preceding section. The only exceptions are the integrals 


é wv 

=f cos fdl = — 

io (re I G3 G 
xf cosafil = [2% — 3% +1 


e€ Ls Sh 4 G G2 +1| : 
=f cos 3fdl = — 4| 2 - OTE | sie Bee 


384 


The separate parts of /:* are as follows: 


constant part: 
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2 og here ee 135 L' 315 L 
Bae 6G 2G Pe. Gr ee 
G? I 10 I 
elec cies. Paes 
sich ee te a 
4 G5 2G® 4, 4 G' 
pa ery ee Bere eh Se es 
es 32°G° 16 G? eer =e 
CO gla Oe aad 
+ Fs 272 aay: ay SH 32 
G 2 
Tp Te ea -: 
soy ls 
eT eh eee ae et 
Wb PP sie s fe ks 
1 SG Paes em aera 
coefficient of cos 2g: 
L$ St b) 
= O15 3 = ee a7 63.L 
4G Ch ens 
a 6 
+ Fi (ore + orn) = +2 — 14 we 
Ree e es: eae 
i ae 4 4 ‘ ACh aa 
LA 
3 Ga PANG yan —2 
ahd ea ieee 
um = 4 eat 
Ap 4° DBs Come 
Co) ta ee pen ie 
G : 2G 2G GL+G 
tse 
3 Ga Pe = of I —2 
a ete edt 
Sum = Che Gl 
coefficient of cos 4g: 
SON jay 0s Es Oi a 
28 64G' ' 32G™ 64G9 
G 5 7 9 
+ 53 272 Gy 32 TK 
G I I 
CT he Saha aie ance 
Sum = 0 
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| The terms with cos 4g disappear, and the expression for /:* becomes, after substitution of the 
| appropriate expressions for A?, AB, B? and BH?/G?, 


f,* 


eee, 2 |) 
Ab BAN: 32 G'® 3 5 G? a2 G's 
a Bk Ee lel 150% gl HEIG 
ee tom) ~ eas gee Ce 


reed A re 2 Hi 
+ Fro | - 3(5 a (: 16 Ga + 15 z=) | cos 2g’, (29) 


in which g has been changed into g’, which is permissible in a part that has k,” as a factor. 


5. Secular and long-period terms. The problem Fot + Fy* ( [SCS ae 1") 


_ is now reduced to the solution of the system of og’? 
| Be orice! equations with the Hamiltonian ee pe pe pee Lope 
2 4 2 
Fe = ar a8 . (- : He : = | + F*, This equation reduces to 
F oy ee F ae 
in which a ce 
Fy* = Fo,* + Fop*, o soa 
i +k * 
_ the former being a function of L’, G’, H only, the en a1 Ue Fee = (33) 
latter depending also on g’. | 0G" og 
Let S* be a new determining function, [byes Wa (34) 
SS = Beye + G''g' = Hh’ Since 
*( TIT a HH US 
“itor ae ig CAG LAE A ys Bi as uke Le A) 
the equation F* = F** may then be written OG 2G" 5 Gn} 


4 in which L’ has been written for L’’, H for H’’, equation (33) gives 


OS\* i ae Ht HH (Ee a ; 
BF ~ rn [5(Be— Be) (o~ B+ 0 2)] sgn) one 


Substituted into 


0S\* 
dg’ ’ 


GH —_ (GY + 


there results 


PAD RAL ek Ht \4 a 
Gane (1 +1 ( Ga ae aa) 3 = Fi an) oa sali a sa | cos 2¢ (35) 


Also 
~~ . We: L"4 I Fe 5 Et H? \- ; m 
Br = G 72\ Ga om) iie\! — en] —2 Gul! — SG ous 
and 
0S\* 
(pei (XD 
ey aL! : 


is 6a Nee IP HS BY] 
= me g(t— gm) —sgali—sgn) |sin2e" - 
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* ree dale aang 55 wee 
=g' +72 16 Gi? I — 33 Gin 56 Gi pe Wek: 
25 L 35 L'*\ Hf m\ 
+(- 22,485, G4 I SG 
L” Ee H® F? —2 its 4 P 
ip oto) G2 aa G4} G's see hs) Gi? sin 2g', (37) é 


0.91" / 
oH es af 


ihe pL 


Le gs 
= w +. (Sa - Ga) 


x [2H + oF ae g(r al 
8 Gr + 10 Gig I 5 Gn 25 Gls aa 5 Gn sin 2g, (38) 


in which g’ has been changed into g’’. 
The Hamiltonian F** is a function of L’, G”’, H only. It is given by equations (32) and (34), 


aie i nee (- Beata 


~ 307! I8G\” 2% 2G? 


) FP DY lege 


where the last term is the first part of the right hand member of equation (29). Let mo be defined by 


Veet 
= 
Lh 
‘Ehen+ 
au _arts 
di ao 
1 C9 Bla 
a {3 Toye aa 2 Fe iz) 
[5 Lee aes gb us POL IEP QUAN V5” 
+ V2 32 Gills 2 G's 32 Coe 16 G5 9 G's i 16G’" J] GG? , 
(5 L227 LE 305 Lee 3 
+ (2 Gis fF 2 G's ze 32 Gi} Gee (39) 
dg" Mt dling > 
dite OGe 
L*( 1.5 © 
meant) {372 an(— 2 ati 2 ie) 
Ws ey Orie TOs Li T8On* dies Neon Oe ) 1BE: 
2 — — —_ ——  _ 
+ 2 E Gi'6 le ae 32 G8 16 G’’é I G7 ate 16 G8) G”? 


13508 135 ES rs Sse 4 
+( 32 GE ls 4 Gua = 32 Gis (Cue ’ (40) t 


ae 


or 
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dh! a F** 


—0~C~*~« MT 


— no | - 
Zing is 
=e y2 |( 2 Gls =F 


TABLE I. FACTORS ARISING IN THE ARGUMENT 
OF THE PERIGEE, EQUATION (36) 


SOKO SIGE 
one G4 a 
oe 
BG) 


Sn) me 
8 G8 Gl 


I 25 cos§ I 5 cost J 
(1 — 5 cos? I)? I —5cos? I 

0° 1.562 —1.250 
10 1.539 —1.222 
20 1.476 —T.142 
30 1.395 —1,023 
40 1.350 —0.890 
50 1.552 —o0.801 
55 2.140 —0.839 
60 6.250 = 1.202 
61 10.58 —1.576 
62 25.72 —2.351 
63 234.76 — 6.956 
64 115.72 +4.716 
65 12.45 +1.120 
66 - 3.790 -+-0.792 
70 0.232 -+0.165 
80 0.001 +0.005 
90 0.000 0.000 


This completes the solution of the problem for 
orbits with inclinations sufficiently far from the 
critical inclination. In the vicinity of the critical 
inclincation the results may become illusory. 

The data given in Table I will serve to show 
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+(- 
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15 L's 
8 (EMG 


a7 L" 
2 (Gun 


105 L’® \ HB 
= 8 Gs] Gs. : (41) 


au" 
dt 


dg” be no| v2(- 
fi ee 

dh" 
COS d= — a no| v2 ( 


vt 
aoe, [r+ = 3 + 12 cos 1”) 


I| 
SS 
o 
ee 
Lael 
= 
S 


3 


By 


=)3 cos? 1’) 


The sidereal mean motion is 


d(l"’ + g’’) x: an” 
Ns = at + cos I di’ 


era) + 


ae 


that for an orbit with moderately small ec- 
centricity, no significant loss of accuracy occurs 
if the inclination is even as close as 1°5 from the 
critical inclination, 63°26’. The second column 
in this table gives the value of the factor that 
contains (1 — 5cos?J)~*? as it appears in g’, 
multiplied by y2e?(1 — e?)-*. In h’ the coefficient 
5 times that in g’. No term of this type 
occurs in l/l’, The third column in the table 
indicates the magnitude of the parts of the 
coefficients that contain (1 — 5cos?/)—“! as a 
factor. In J’ and g’ these numbers are multiplied 
by +72 times a factor that reduces to unity for 
é = 0;inG’ and hf’ the coefficients are diminished 
by a function of e that contains e? as a factor. 


6. Comparison with results obtained by Atll’s 
method. For e = 0, the motions obtained in 
equations (39), (40), (41) may be compared with 
my results obtained by Hill’s method (Brouwer 
1958). With L’/G”’ = 1 these expressions reduce 
to 


+6cos I” — oF cost |, 


a, —~ — 36cos? I” + oe cos! 1’) | (42) 


I 28 
Ns = No E eye 3-1 9 cos”_T'’) “+ v8(+ r SoRGOLCOS cl casts = cost r’)|, 


I 
io Me [: + yo(+ 3 — 9cos? J”) + y2? (+ — 24 cos? I” + cos! ”)|, 


mn —" 


F] 
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i 
Thus | 
7h Ee My a 
za = gs Ns E +b al -F 3 cost) rata (+ S.Cos! ae 5 cost |, , 
i 
A dh" / ‘yt 3 vt | 

cos I Thee v2(— 3 cos? I’) + 722 { — 3.cos? I — 5 cost I ; (43) 
In the Hill method applied to the artificial Now a 
satellite problem I found (Brouwer 1958, p. 435) Gy ante j 


d(l!t + ¢’’) 
dt 


= NsH E ap 372 cos” Ir 


+ y,' (2 cos? Ix — 30 cos! I) | (44) 


the subscript H being used with ms, J. 

The expression in square brackets is equivalent 
to the value of 1/(1 + p) for ys = 0 with 72 
replaced by vy’. In order to compare with the 
result obtained above it should be noted that 
yo', cos Ig and msy differ from y2, H/G"’ and ng. 

In the present theory, if L’ = (uao)?, 


wko ke ger | spt 


jy = Sar nH =. — 
q ee ao : i? ae 
In the Hill theory 
ities ke re a 
s \) Ss a? 


Hence 


2 
w 
I| 
Q 
Se <) 
ll 
Fa 
Sa 
feo] co 
~~ 


l| 


I + 472(1 — 3 cos’ In), 


from equations (42) and (44). Hence 


(45) 


The inclination Iq in the Hill‘method was 
defined as the maximum latitude. In the present 
theory this does not apply to the inclination 
obtained from cos /’’ = H/G” in view of the 
presence of the term of the first order in y. and 
zero order in e in G/G’. For e = 0 the relation is 


G = G[1 + 37 aed cos (2g + 21) ]. 


yo = ¥2 — 4y2(1 — 3 cos? Iq). 


Then 
PIECING 
Acos? J = — G2 Gr 
= — 3y2 sin? I cos? I cos (2g + 21) 
or 


sin? J = sin? I’’[1 + 372 cos? I” cos (2g + 21) ]. 


d(l!’ a g'’) 


sino = 
‘= gin? J’’[1 + 372 cos? I” cos (22 + 21) ae 

x sin? (g + 1). | 

This isa maximum for 2g + 21 = 180°, when | 
sin? Bmax = sin? Jy = sin? I”’[1 — 372 cos?’ aa 

from which ; 
cos? In = cos? I/"[1 + 3y2sin?I’’]. (46) 


The introduction of (45) and (46) into (44) © 
yields 


t 


= NaH [: + 372 cos? I” 


Introduction into (47) yields 


oa 


dt q 

9 15 92 vt 4 vt ; | i 

ae Ye (= cos: Eo acoso Ste i 

= i 

Finally * 
d(I!’ ae g’’) dh!’ | f 

1a = ai + cos I ae f 

d(l!’ me g'’) dh" 4 

NsH *= it + cos In di’ i 

or a 
wt id 

Nsw — ns = (cos In — cos I’) a ¢ 

a Ww 2 

ee sin? W/icosue” a 

=— : nsy2 sin? I'’ cos? I”. j 


d EG aA . 
mee es Aly + 372 cos? ype” 


+ y2?(3 cos? I” + $ cos? oe 


in agreement with the expression (43) obtained — 
previously. 

In addition to establishing the agreement 
between the results for e = 0 obtained by two 
different methods, the discussion serves as an 
illustration of the caution necessary in the com- 
parison of results to the second power of v2 
obtained by different treatments of the same 
problem. 
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7. The fourth harmonic. The coefficients of these terms are so small that it is unnecessary to consider 
the terms of short period. Hence the contributions to the Hamiltonian may be limited to the parts 
independent of /, and only the first powers are needed. 

The expression for U, is 


wks 


Ty = Ab — 10sin? 6 + ia, a). 


y? 


Since 


I 
Sipe. 5 sin’ I jry— cos(2¢ ---2:f)4, 


Syl (Sh = sin! I([3 — 4cos (2g + 2f) + cos (4g + 4f) ], 


this becomes 


lL 8 4 


5k 5 5 ? 20 a? 
U,=" £1] (3 = S cost 1 + Seost 1) 5 + ( ah z cos? I 35 cost 1) © cos (2g + 2f) 


+ (= rp £08 Hee cos r) 75 COS (4g + 4f) }. 


The secular part is 


wR See 635 Ht 5_ 33) 
LIGTINS AG oy (Ge Demerow 


The part independent of g is added to F;,*, 


haf isl’ 9. L ( HP 35 eS 
mee ee ee oety\ OG? | 3 G4y? 


5 20H? sZ)(3 5S) co | 
S 2 : 
Ea Br 


aE 

laa 
| 

nie 

o® 

fon 

Q 


AaFs5* 


the part having cos 2g as a factor is added to Fop*, 


: ral Get Ee ee IEE ¢ 
AaFop* = : a Gi Ga JE 8 ek. Sai Gia ) COS 28. 
Now 
0 QeleG 4 Jiebo NE 
eve = 22 (1 = Sa) AgFoy* 


oe Gh es eee te ee 
= 5 Ea" (Za - Sa) (1-8 Bat 7G) (1-52) cos 2g 


ce ke, (Le 1? iP HT ee 
Bee (an - ae I — 3 Gm — 8 Gm ee ON eIp cos 2¢’. 
By integration, 


2 k Ie Ee FH? Hf H? tl : : 
AgSi* = 5 a = Gr (4a ba =) E TE ie) G'? = 8 Gi'4 (: 5 ay | SU eee 


If now 
utks 
i L’8? 

so that 
wks DORE 
Lk» uy 2 


aoe al 
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it follows that the additions to the right-hand members of equations (35) to (38) are 
a. a ple a fe ALE og Dg ea ee z) | if 
AyG’ = eine G (5 awe I 3 Gi? 8 Gi I 5 Gin cos 2g”, 
1) Lares HT” an r aa) | Py 
Ad) See 12 Yo G" [: 3 Gi? 8 Gls I 5 G’”? sin 2g , 
7 ee sz 5 An) - & al di a) (- 85, 25 Z) 
Aug = v9 E G4 I 5 G2 24 Gi? I 9 G’? =F 3 Gi4 = 3 Gi’? 
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The additions to the secular motions are 


| a eee a i( 2 Ge ue Ne see 
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27 Are G6 — 35 Eis G4 Dy, 
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8. The third and fifth harmonics. Odd harmonics are to be considered only if symmetry of the earth | 
with respect to the equator does not exist. A study of the motion of 1958 62 led to the introduction 
of these odd harmonics by O’Keefe, Eckels and Squires (1 959). 


The expression for U3 is 
HA s.o tae nilt ome 
U3 = ee (: sin’ B > sin a). 


With i 


sin? 8 = sin? I | Sin (g+ f) — zsin (24h 3/| 


this becomes 


‘ 


? : 
bee ale sin r+ ® sin’ 1) sia (gis) a 2 sin? Tsin (3g + 3/) | 


v4 


For the sake of brevity, let e’ sin J’ stand for (1 — G”/L”)#(1 — H?/G”)}, e’’ sin I for the same 


expression with G’ replaced by G’’. The secular part may then be written j 


5A key 
A3Fop* = aa e’ sin I’ (3 — 3) sin g’. * 


There is no part independent of g’. 
It follows that 


) I A3.9 L” 
ae! 1% = Tune 7, Gh Gin Ge" sin I” sin g’; 
I Aso L”? 
At = -- 5 = = Ge sin I’ cos 9". 
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i if 


With 


| the following results are easily obtained : 


Bi at INAH 


/ 


A3g => 


A;3h’ = 


The fifth harmonic is 


Us 
With 


this becomes 


uA 5.0 iS ore 
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Var L's ’ 
I LP 
A;G’ = —-— AG Gin Ge” sin I’ sin g’’. 
; sin I G? 
aL @sin i) = ae 
my au sin I G G@ Jee 
CS as ae e D2 sinI G 
fo) : G tél 
Pee inl @ 
I 3 sin EGS 
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472G e sin [’’G 
ile IGS J66 Qe i 
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Syl (6) sine 1 S sin (¢ +f) - > sin (e + 3f) + = sin (5g + sf | 


(0 Nes ae sin? 8 +3 2 sin’ a|. 
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| 7A I Ir TsO, 
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| The secular part is 
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) Sain (el eiaesif ho °3. sin’ 1S sin (5g + sf}. 


)z sin (5¢ + sf) |. 


G 
) ( — =) sin 3¢' | 


— 


' 
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Hence 
Ce) , WA 5.0 Je : 4 
at ee gus eae egret ’ 3 
ag AS GisR, (: 5 a) CL Cina F 
ue EP iH )( + =) ; j ; 
x | - 64 ( 14 7m G'? 21 Gi’ i 3 L” sin g 
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the resulting long-period terms become: 
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37 3 + 16 Gi? by 


15 
‘: [38(7- 


35 ( 
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9. Formulas for computation. For convenience of 
computation the perturbations in the Keplerian 
elements a, e, J are given instead of those in 
Gall. 
The adopted force function is 


es asi 
= tg (WF — 3sin’ 8) 


‘ 


k 
+ (1 — 1osin? 6 + 35 
r 3 


A 
ak ee (-3 3 sin 8-B 2 sint s) 
A 
+ ease sin 8 eae + Sains), 


in which kz is a small quantity, and ka, A3.0, As.o 
are assumed to be of order k,?. 
The secular motions have been computed to 
O(k.”), the coefficients of periodic terms to O(k»). 
Basic constants: 


a’ = semi-major axis constant 
e’’ = eccentricity constant 
I’ = inclination constant 
mo = wa" 4 = 17.04337 (do/R) rev./day 
R = equatorial radius 
Abbreviations: 
m=(1—e'”)t #6=cos I” 
Re Ra A 3.0 A 5.0 
Ve ia va Bahai la 5 q's 
v2 =o ya =yan® vs’ = an 8-8’ = Yon 


Secular terms: 


I’ = “mean’’ mean anomaly 


= not 


nee 
ae 5 32 G2 I 


“1 Ht a 
san | + 40 Gals - sz) | cos @” 
Er es (CS EY" ccs ae” 
5 Gi? e) G4 I 5 Gi? COSNS See fe 


It is customary to use for the second harmonic 
the coefficient J; Jeffreys (1954) used for the 
fourth harmonic the coefficient D. The relations 
between J, D and yo, y4 are 


I ig We 2 Ike Ne 
oes wt ) ee al! : 


Strictly speaking, e’’ + 6,e, 6’ = cos I” + 6,1), 

= [1 — (e’’ + 6,e)? * should be used in the 
computation of the periodic terms, but since the 
short-period terms are obtained to O(k:), it is of 
no consequence if contributions of O(R»2) are 
omitted in expressions that have 2 as a factor. 
Similarly, 1’, g’’ might be used in computing f’, 
r’; but since l’, g’ are available, their use does not 
complicate the calculation. 

The formulas are applicable for any eccentricity 
e < I and any inclination with the exception of 
inclinations near the critical inclination, for 
which 1 — 5 cos? J appears as a small divisor. 

The appearance of e’’ as a divisor in the short- 
period terms in e is apparent only. The expres- 
sions that are multiplied by e’’ contain e’’ as 
a factor, either implicitly or explicitly. 

In the short-period terms in / and g a divisor 
e’’ occurs also, but for the calculation of the 
position only g +/+ equation of the center is 
needed. In g + / the divisor e”’ is not present. 

Singularities in some of the elements also 
occur for very small inclinations; again, no 
singularity is present in the coordinates. In such 
cases it may be found convenient to modify the 
formulas and obtain expressions for the per- 
turbations in coordinates. 


i 3 y!n(— a = y2/*mL— 15 + 169 + 250? + (30 — 960 — 907”)? 


I 
+ (105 + 1440 + 250°)6*] + = ya'ne’"[3 — 308 + 356°] ¢ + Lo’ 
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g’’ = mean argument of perigee 


= not | S y)(— 1 + 50) + = v2"L— 35 + 240 + 250? + (90 — 1929 — 12672) 


+ (385 + 3600 + 4517)64] + Fe > y'[21 — on? + (— 270 + 1267n?)6+ (385 — 1891904] + go” 


h'’ = mean longitude of ascending node 


I| 


not ! = $700 $ "E(- 5b 129 -997)0 4 = 35) eg = Bae | 
“i = yl (5 = Bar CS see 70)| + ho” 
Long-period terms: 
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I 
W=h" + {- g vee Lit + 8007(1 — 56")-! + 20064(1 — 562)-? ] 
by v4" 
I2 Y2 


Mya es @ 5 v5 e6 
472 sinl” § 64 72’ sin I” 4 


Ya ingr + 1667(1 — 567)! + 4064(1 — 562)-?]; sin 2g” 


=eiGen a) EIU 25240 (1A 56") 2] 


K5 15. 


Stead aa e”6 sin I” (4 + 3e’”)[3 + 1662(1 — 567) + 406411 — se] cos g”” 
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Short-period terms included: 
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cos (2g + 21) |} 
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tye 0 ize cos (2g + f’) + e” cos (2g’ + 3/1] 
T= I" + 6,0 + 4y2'0(1 — 6°)3[3 cos (2g’ + 2f’) + 3e’’ cos (22’ + f’) + e” cos (2g’ + 3f’) J 


. ni . ee al’ : 
ep ree one: 2(— 1 + 36?) e+ = 2) sin Ff 
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+3¢1 — 0 |(- a Be}. (2g +n+(¢ me +5 eee ~) sin (2g +3/]| 
ene tan fata (Get trains 
4e” 2 7 1 y! 


+ 3(1 — &) I(- ao :) sin (2g + f’) + (= ae =) sin (2¢’ + 3/)|| 
r r r r 3 
+ ty/{6(— 1 + 50?)(f’ — + e” sin f’) 
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h =h' — 3y,/6[6(f’ —’l’ + e’ sin f’) — 3 sin (2g’ + 27’) 
— 3e’’sin (2¢’ + f’) — e” sin (2g’ + 37’) ]. 


f', r’ are to be computed from 
El —e"sin EH’ =I 


Tite: 3 \2 is: 
tan 3f’ = i tan 3F’ aren f = (i — e’)* sin E’ 


Mt iad 7 / 
a 1 + e”’ cos f if 
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For the calculation of the coordinates at any 
time the complete values of e and / should be 
used for the solution of Kepler’s equation, 


E-—esink =] 


and subsequently r and f, which may then be 
used in the formulas: 


r[cos (g + f) cosh — sin (g + f) sink cos I] 


\| 


x 


y = rLcos (g + f) sink + sin (g + f) coshcos I] 
z=rsin (g + f)sinI 
A convenient alternative form is: 
x = Ax (cos # — e) + B, sin E 
y = Ay (cos E — e) + By sin E 
z=A,(cosE —e)+ B,sinE 
Ax = a[cos gcos hk — sin gsinhcos J | 
B, = — a(t — e?)F [sin g cosh 
+ cos g sinh cos I | 
Ay = a[sin gcoshcosI + cos g sinh ] 
By = a(1 — e)# [cos g cosh cos I — sing sinh ] 
Laplace Be Bs 
H. Struve = ze = 
W. de Sitter Came 4 KR 
a 15 
D. Brouwer — 2ky sh 
Hovetrers = 2K 3 pp 
35 
Y. Kozai eerie ays 
35 
P. Herget and P. Musen — 2k» 8 ka 
J. O'Keefe et al. + Azo/u + Aso/p 
B. Garfinkel — 2k k’ 
J. Vinti — JR? — J4R! 


In the table RX represents the earth’s equatorial 
radius. Ignoring the presence of R? and differ- 
ences in sign, essentially three different coeffici- 
ents for the second harmonic have been used in 
recent papers. For the coefficients of the fourth 
harmonic six different choices are listed. I now 
regret that I introduced ks, k4 in my paper in 
1946. The principal reason was that they give a 
particularly simple form for the expression of the 
potential in the equatorial plane. If I could have 
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AA 
B, 


asin g sin I 


I 


a(1 — e)} cos g sin I 


Noted added in proof. The lack of uniformity in — 
notation of the coefficients of the second and — 


fourth harmonics of the earth’s potential in 
papers dealing with the motion of artificial satel- 
lites calls for a comment on this subject. 

The table below contains a listing of some of 
the designations used and their relations to the 
coefficients B, in the expression of the force 
function of a body with rotational symmetry, 


pe aE me > oe B) i} 


p=2 


in which P, are Legendre polynomials and yp 


= GM. The expression is an adaptation of the 
Laplacian expression given by Tisserand. 

In addition to the equivalents of B, and B, the 
table gives those of the ratio B,/B,”, which is 
unity for the special case treated by Vinti (1959), 
in which the terms with small divisors near the 
critical inclination vanish. No effort has been 
made to make the tabulation complete. 


Bs/B? Tisserand, Méc. Cél. Il, 
320, 1890 
3 1/k Suppl. I, Obs. Pulkovo 1888 
oK/S: B.A Ns2, 07 p19248 
5 balk! A. J. 51, 223, 1946 
18 
peie M. N. 14, 433, 1954 
18 
mea et 
ag 
2ka/ko? A. J. 63, 430, 1958 
pA 4 o/A2, o? A. J. 64, 235, 1959 
ele This issue 
Jif FP J. of Res. Nat. Bureau 


of Standards 62B, 105, 1959 
f 


foreseen the increase in interest in.the subject 
and the confusion to which I was contributing, I 
would have chosen the coefficients B, or the 
alternative form 


r=*[1-53,(4)'2, cno)], 


which was used by Vinti (1959). I intend to revert 
to this form and recommend this to other authors. 
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Anyone who has ever carried out develop- 
ments of the type presented here knows that 
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extreme caution is necessary if errors are to be 
avoided, even if the method used is straight- 
forward. The comparison with Dr. Garfinkel’s 
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PHOTOGRAPHIC DIMENSIONS OF THE BRIGHTER GALAXIES 
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Abstract. Standard photometric dimensions of 37 galaxies are derived from recent detailed photometric studies. The 
relations of photometric dimensions to the laws of luminosity distribution in galaxies of different types are discussed. 

Micrometric dimensions of over 1200 Shapley-Ames galaxies (to be published later) have been derived from 13 main 
series of photographic determinations, including 6 series of new data. The relations between systems of micrometric 
diameters D and axial ratios R = d/D are determined and the relations between photometric and micrometric dimen- 
sions are discussed. 

A provisional reduction of 10 series to the (corrected) system of the Heidelberg survey by Reinmuth was made for 
the Shapley-Ames objects. The probable errors of the reduced values from individual series are of the order of 12 per cent 
for both D and R. 

Strong correlations between apparent diameter and axial ratio, i.e., tilt, and galactic latitude are indicated. The re- 
duced or ‘‘face-on”’ diameters D(o) of galaxies of a given type are also closely correlated with total magnitude mr and 
can be used as distance indicators. An application is made to the calibration in terms of the current distance scale of the 
relative distances of several clouds and clusters of bright southern galaxies observed in the Mount Stromlo survey. 

The application of the corrected mean diameters and revised luminosity laws to the determination of improved incom- 


pleteness corrections and total (or asymptotic) magnitudes is briefly discussed. 


In the preparation of a reference catalogue of 
bright galaxies based on the Harvard ‘‘Survey of 
Galaxies brighter than the 13th magnitude’”’ 
(Shapley and Ames 1932, de Vaucouleurs 1953b) 
auxiliary studies have been made of the magni- 
tude and dimension systems used in various 
published and unpublished surveys in both hemi- 
spheres. The discussion of photographic magni- 

tudes published three years ago (de Vaucouleurs 
1956a, 1957a) has indicated the need for a de- 
tailed study of dimension data. The present paper 
summarizes the main results of this study and 


‘ 


discusses briefly some applications; further de- 
tails will be found in Annales de l’Observatoire du 
Houga, Vol. II, Part 2. The revised reference 
catalogue is now being prepared for publication. 

1. Photometric dimensions. The apparent di- 
mensions of galaxies can be precisely defined and 
correctly determined through a detailed photom- 
etry of the luminosity distribution in the nebular 
image. The brightness dimensions D(m,), d(m1) 
refer to a given isophote (or best fitting ellipse) 
of specified surface brightness, for instance m, = 
25.0 mag/sec? (pg) (Redman 1936; Hazen 1957). 
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The effective dimensions D,., d, refer to the 
isophote (or best fitting ellipse) within which is 
emitted half the total (photographic) luminosity 
of the nebula (de Vaucouleurs 1948, 1953a). 

The determination of brightness dimensions 
requires absolute measurements of surface bright- 
ness on a given magnitude scale in mag/sec’; 
these dimensions vary with the particular bright- 
ness level arbitrarily selected in their definition. 
The determination of effective dimensions re- 
quires only relative measurements of surface 
intensities with arbitrary zero point; the dimen- 
sions are, in principle, intrinsic properties of the 
nebula, independent of the observational circum- 
stances (except, of course, wave length of light). 

Both types of measurements are slow and 
tedious so that reliable data on either brightness 
or effective dimensions are at present available 
for only a few dozen objects listed in Table I. 
More work along these lines is greatly needed. 

On the other hand, ‘‘maximum’”’ dimensions 
Dax, max Measured on microphotometer trac- 
ings are now available for several hundred bright 
galaxies observed under fairly homogeneous con- 
ditions; these values are approximate brightness 
dimensions relative to an unspecified mean (pho- 
tographic) brightness threshold, usually fainter 
than 25 mag/sec?. Comparison with calibrated 
photometric data determined this threshold to be 
about 25.3 mag/sec? (Patterson 1941) and 26.7 
mag/sec? (Holmberg 1958). 

All photometric dimensions which are in the 
nature of brightness dimensions are determined 
by the threshold brightness B; and the law of 
luminosity distribution B(r) in the image of the 
galaxy, while the effective dimensions depend 
only on the latter. The ratio of brightness dimen- 
sions relative to different magnitude thresholds 
My, M2, --: depends on the luminosity laws; ex- 
perience shows that, in practice, only two main 
groups need be considered in a discussion of 
apparent diameters: 1) the ellipticals, 2) the 
lenticulars and spirals. : 

The first group obeys closely the luminosity 
law log B = ar/4+ 0 (dé Vaucouleurs 1948, 
1953a, 1956a), which may be written in the re- 
duced form 


log @ = > 2. 33(os 1) (1) 


if @ = B/B., a = r/re, te, Be being the effective 
radius and corresponding surface brightness. 
The second group follows approximately, ex- 
cept near the center, an exponential law, log B = 
A — Kr (Patterson 1940, de Vaucouleurs 1958), 
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TABLE I. PHOTOMETRIC DIMENSIONS OF 37 GALAXIES 


NGC Type De de[De Di* di/D1* Sources 
LMC SB(s)m 5°5 0.89 == = I 
SMC_ SB(s)mp PScn On5% = —= I 
224 SA(s)b D2 OnSO: 3°05 0.29 2 
598 SA(s)c 24S) (OCONEE eue 0.62 g 
1291 SB(s)O* 7B Ow 5S aad re 4 
1313 SB(s)d 4600) JORG. = = 5 
1316 SOp Silene =2e ae 5, 6 
STS) en 3(p | 10.32) a a4eOns 7,8 
3379 EO Ps | ARO = = 8.4 
4267 SB(s?)O7 — = Pa 0.96 9 
4350 . SAO sp 1/5 0.50 B02 0.62 9 
4365. E3 PIR SONS O47) EROnyh 9: 
ABT Ae Et 3/0 . 0.90 Fon, 0.99 9 
4406 Et, 646. 0170. (2.2) aonae 9 
4417 SO 146 0.45 315 < 2ORA0 9 
4425 SO sp. — — Ay 0.40 9: 
4442 SB(s)O° — == 4:6 0.40 9: 
4459 SA(r)OF ae a 4:2 0.83 9 
4461 SB(s)O?: — == 3:5 0.40 9 
4472 2 — — 1276)" SOnSO 9 
4473 1/9 0.60 417 0.63 9: 
4477) USB(st)O2;) — 2:6) 20200 9 
4486 Eip 2'8 0.88 gil 0.75) 16559) 
4494 Er Te2 ue OLO': 5/5 COMO Mayes. 
4503 SBO- ae ese 3:5 0.47 9 
14520 “SAB(s:)O> — — Tas O37, 9 
4552 EO 129) On90. 5/6 0.87 9 
4564 SO? 1‘6 0.50 3.0.) Ouse 9: 
4570 SO sp a ni 3:8 0.33 ce) 
Mephele oO): — — ai 0.70 9 
4504 SA(s:)a 4462210250 = — 10 
4621 \ E5 3.0), 10570 5/7 0.82 9 
4649 E2 250) FORSO 9/2 0.78 9 
4694 SBOp = a 3:4 0.57 9 
Topsils ) pele ngcey (OF 9 = =F 4:6 \Oags 9 
5643 SA(s:)c 23 aD 300 — — 5 
6744 SAB(r)be OR 0508 == ae 5 
*m, = 25.0 mag/sec? (pg). 
Sources: 
1 = de Vaucouleurs 1957b 
2 = de Vaucouleurs 1958 
3 = de Vaucouleurs 1959 
4 = de Vaucouleurs 1956c 
5 = unpublished 
6 = Sersic 1958 
7 = Oort 1940 
8 = de Vaucouleurs 1953a 
9 = Hazen 1957 
10 = de Vaucouleurs 1948 


which similarly can be written in the reduced 
form 


log ® = — 0.729 (a — TI) (2) 


with @ and a defined as above. __- 

There is @ priorz no constant relation between 
brightness dimensions and effective dimensions 
since it depends on B, and luminosity law; how- 
ever, if the dwarf ellipticals and irregulars are 
excluded, the mean surface brightness of the 
majority of catalogued galaxies is restricted to a 
fairly small range (cf. de Vaucouleurs 1957a, 
Fig. 19 and sections 5 and 6 below) so that, for 
a given type, a statistical relation exists between 
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these dimensions. For the small sample in Table I 
brightness dimensions relative tom; = 25.0 mag/ 
sec? are, on the average, about 2.5 times the 
effective dimensions: D(m)/D.'~ 2.5. It is then 
possible to estimate roughly, by means of equa- 
tions (1) and (2), the variation of brightness 
diameter corresponding to a given variation Am, 
of the threshold brightness when m, ~ 25; for 
Am, = +0.25 mag, Aloga = +0.04, or IO per 
cent on D,, for ellipticals, and Aa = +0.14, or 
5 per cent on Dy, for spirals. The exponential 
luminosity law in the outer parts of spirals makes 
their diameters less sensitive to fluctuations in 
the brightness threshold; this is confirmed by a 
direct comparison of micrometric diameters (sec- 
tion 4). 

2. Micrometric dimensions. For the great ma- 
jority of galaxies the only dimension data avail- 
able are direct micrometer measurements of pho- 
tographic negatives. These micrometric dimen- 
sions refer to an unspecified mean brightness 
threshold depending in an a@ priori unknown 
_ manner on the properties of the object (size, 
shape, luminosity gradient), of the foreground 
(airglow, galactic and zodiacal light), of the 
photographic plate (exposure, graininess, con- 
trast factor) and of the observer and observing 
conditions. However, when the plates are secured 
and measured under homogeneous conditions, 
micrometric dimensions refer to a relatively well- 
defined brightness threshold which can be cali- 
brated by appropriate comparisons with photo- 
metric data. 

The main purpose of the present paper is to 
investigate the systematic errors to which micro- 
metric dimensions are subject, to obtain the re- 
duction formulae to a standard mean system 
and to calibrate the mean system of brightness 
dimensions so defined... 

The material available for a study of micro- 
metric dimensions of bright galaxies is summar- 
ized in Table II. Altogether over 30 lists and 
catalogues grouped in 13 main series were used ; 
series IV, V, VI, VII, VIII, LX include new mate- 
rial obtained from the writer’s examination of the 
following plate collections: (IV) Isaac Roberts 
20-inch, now at the Paris Observatory, (V, Va) 
Mount’ Stromlo 30-inch Reynolds Survey (de 
Vaucouleurs 1956b) and 74-inch survey (de 
Vaucouleurs 1960), (VI) Lick 36-inch Crossley 
(modern plates only), (VII) and (VIII) Mount 
Wilson 60-inch and too-inch, (1X) Palomar 200- 
inch. In series (Va), (VII), (VIII) and (IX) two 
sets of dimensions are available for each galaxy: 
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D;, d; for the inner bright regions and D,, d, for 
the outer, fainter regions (cf. de Vaucouleurs 
1956b). 

After some experiment it was decided that the 
significant parameters for a comparison of differ- 
ent lists were the major diameter D and the axial 
ratio R = d/D (not D and d separately ; see also 
Holmberg 1946). Further, the discussion was 
limited to objects listed in the Shapley-Ames 
catalogue, which offers the greatest overlap be- 
tween the ‘various series. 

The comparison of the dimension systems D,, 
D: of two catalogues was made graphically on 
logarithmic correlograms (Fig. 1 right), and be- 
cause it includes the largest number of objects 
north of declination — 20°, the Heidelberg survey 
(series Ia, Reinmuth 1926) was taken as the 
standard for comparison purposes. In general an 
equation of the form 


D» = AD, (3) 


suffices to represent the relation between series | 
and any other series (II, III, IV, V, VI, etc.) of 
dimensions secured with small telescopes. How- 
ever, when the comparison is with dimensions 
given by the larger telescopes (series VII, VIII, 
IX etc.) an equation of the form 


D,=AD,+B (4) 


is often required. The reason for this difference 
will be given in section 3. Experience shows that 
ellipticals on the one hand, lenticulars and spirals 
on the other, give different coefficients in (3) or 
(4) as could be expected from their different 
luminosity laws (section 1). In a few cases a non- 
linear relation was indicated, pointing to a faulty 
dimensions system in some series, especially for 
ellipticals (series IIb, XIb). 

The comparison of axial ratios was made simi- 
larly on logarithmic plots (Figure 1 left). In 
general, two groups of series are indicated: (a) 
series for which the axial ratios R are nearly on 
the same system as that of the Heidelberg survey 
(series Ia) whose systematic errors have been 
determined by Holmberg (1946); and (b) series 
for which the axial ratios are in agreement with 
those determined by photometry and whose sys- 
tematic errors are consequently small or nil. 
This grouping was found by comparing the em- 
pirical (Ri, Re) relation to the relation predicted 
by Holmberg’s corrections. A relation distinctly 
different from that found in either (a) or (b) was 
observed in only a few cases (series VIII, XII, 
XIIIb). 
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TABLE II. SOURCES OF MICROMETRIC DIMENSIONS OF GALAXIES 
Observatory and Scale Number of galaxies 
Series telescope Imm = Total In HA88, 2 Sources 
la Heidelberg, 16’’ a7 3500 776 I 
Ib = a= 827 149 2 
Ic == — 4866 426 2 
Ila Harvard-Boyden, 24” 1/0 243 243 4 
IIb ns = == 486 486 5 
IIla Helwan, 30” 0/9 58 40 6 
IIIb — = 184 QI 7 
IIIc — — 169 60 8,9 
IIId a= = 461 65 10 
IV Roberts, 20’ Tees 600 _ 289 II 
V Mt. Stromlo, 30’ 1/12 460 205 12 
Va Mt. Stromlo, 74” 0/38 47 40 13 
Vb Pretoria, 74” 0/38 34 19 II 
Ve Harvard-Boyden, 60’ 0144 IOI 61 II 
Via Lick, 36” 0/63 494 105 II 
Vib al == 54 Il 
VIc = a= 59 1 
Vila Mt. Wilson, 60’ 0145 384 68 II 
VIIb = = 166 Il 
Vile = == 53 II 
Villa Mt. Wilson, 100’” 0127 503 51 II 
VIIIb = —= 182 II 
VIlIc == == 167 II 
IXa Mt. Palomar, 200’ 0/20 262 118 II 
IXb = aa 57 Il 
x Lick, 36” 0/63 477 257 14 
Xla Mt. Wilson, 60’ 0'45 IOI 29 15 
XIb Micke 367 0/63 400 303 16 
Mt. Wilson, 60’’, 100’ 0145, 0/27 
XII Franklin-Adams, 10” 3/0 122 96 17 
ick 36" 0/63 
XIa Helwan, 30” 0'9 202 162 18 
XIIIb Lick, 36” 0/63 
XIIIc Mt. Wilson, 60’’, 100’” 0/45, 0/27 
Sources: I = Reinmuth 1926 : 2 = Holmberg 1937 
3 = Reiz 1941 : X(¢ 4 = Shapley-Ames 1932, ‘‘authority (a)’’ only 
5 = Shapley 1934 6 = Knox-Shaw 1912, 1915 
7 = Gregory 1920, 1921 8 = Know-Shaw 1924 
9g = Madwar 1935 10 = Svenonius 1938 
11 = de Vaucouleurs unpublished 12 = de Vaucouleurs 1956b 
13 = de Vaucouleurs 1960 14 = Curtis 1918 
15 = Pease 1917, 1920 16 = Hubble 1926 
Lyfe 18 = Danver 1942 
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The coefficients A and B of equations (3) and 
(4) and the axial ratio group (a or b) are listed 
in Table III; both the raw empirical values and 
the adjusted values resulting from the discussion 
in section 3 are given. The scale coefficient A 
(for the series taken with the smaller telescopes) 
is correlated with the mean epoch of the plates; 
prior to 1910-20 A ~ 0.9 (E) and 0.95 (S) as 
against A ~ 1.3 (E) and 1.15 (S) since 1920-30; 


this is clearly the effect and measure of increased 
plate sensitivity and fuller exposure. 

The sum A + B, measuring the diameter D» 
which corresponds to D(Ja) = 1/0, is correlated 
with the plate scale (series earlier than 1920 are 
excluded) ; a galaxy of diameter 1‘0, as measured 
on small scale plates (0/5 to 1!0/mm) is 1/5 on 
medium-scale plates (2’ to 3’//mm) and 2’ or 
more on large-scale plates (4’/mm and over). 
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This variation is brought about by the relation 
between linear image size and threshold bright- 
ness or, in other words, by the signal/noise ratio 
associated with the graininess of the photographic 
image, as discussed in the following section. 

3. Relation between photometric and micrometric 
dimensions. The threshold brightness is essen- 
tially determined by the smallest perceptible 
excess of photographic density above background 
that can be significantly associated with the neb- 
ular image; this significance is judged by com- 
parison with the density fluctuations of the plate 
background determined by the graininess of the 
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photographic image. Since the relative density 
fluctuations increase as the scanned area de- 
creases, the density excess that can be considered 
as significant increases also and in consequence 
the threshold brightness must vary with image 
size. This is confirmed by observation; for ex- 
ample, Figure 2 shows the variation of threshold 
surface brightness as a function of the diameter 
of uniform circular images as empirically deter- 
mined by Hubble (1932). The threshold bright- 
ness is independent of image size only when the 
diameter is I cm or greater; for a I mm image it 
is about 1 mag brighter. 


ve Ta 


Figure 1. Comparisons of micrometric dimension systems: Ia (Reinmuth 1926) and Ila (Shapley and Ames 1932) 
Correlation between diameters D (right): the scale ratio D(Ia)/D(Ila) is given for E, SO, and S objects. Correlation 
between axial ratios R = d/D (left): the curve represents the relation according to Holmberg (1946) between Rein- 


muth’s data and the true (photometric) axial ratio. 
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TABLE III. RELATIONS BETWEEN MICROMETRIC DIMENSION SYSTEMS 
Ellipticals: A + B Spirals: A+B Axial ratio 
Series Empirical Adjusted Empirical Adjusted E Ss Notes 
la I .00 1.00 1.00 1.00 Ia Ia Std. 
Ib 1.00 1.00 1.00 1.00 la la 
Ic 0.80 0.85 — 0/05 0.90 .92 — 0105 = * (1) 
Ila 1.18 0.95 + 0/2 1.07 0.9 +0/15 ptm ptm 
IIb 2 — 1.43 1.25 + 0/35 ptm. ptm (2) 
III — — 10; 0.85 + 0/15 — la 
1V 1.05 0.95 + 0/1 1.00 0.95 + O'1 Ia la 
VD; 0.94 0.85 +0!1 0.84 0.80 + 0/05 ptm ptm (3) 
—Dp 1.40 [20 = 023 1.20 I.I +0(2- ptm ptm 
Va D; 1-56: 0.8 +043 0.83 0.65 + 0/2 ptm.? ptm 
— Ds 1.5: +1!0 1.4 +0!6 esa 1.2 +05 ptm.? ptm 
VI 107: 0.9: +0/3 1.06 0.85 + 0/2 la? la 
VII D; 0.84 0.6 +0!2 0.82 0.6 +0/2 la la 
—Dy 1.4: +0%3 1.2 +0!5 TOs aOn st ti a3) oA) la la 
VIII D; 0.6: 0.6 +0!/2 0.82 0.6 +02 ptm.? ~~ la 
== Bhp 1.4: +04 1.3 +0/6 1-10 +07 1.15 + 0/5 ptm fa 
IX D; 0.8 +0!2 0.7 +0%3 0.8 +0/3 0.8 +035 ptm. ilar 
—Do 1.8 + 1/0 1.9 +0!9 [3 re, 1.4 1:07 Ia? la 
x 0.77 — 0.95 = la Ia 
Xia 0.9: — TO — la Tat : 
XIb és — 0.95 — ptm la (4) 
XII 0.82 - 0.95 = : Ia (5) 
XIIla — — 1.15 = a la? 
XIIIb -- — 1.10 a an _ (6) 
XIIIc =e = 1.00 = = Ia 
Notes: 
(1) after Holmberg’s discussion (1946) R(Ic) = R(Ia) — 0.02. 
(2) the D scale for E is non-linear: log D(la) = 0.86 log D(IIb) — 0.17. 
(3) coefficients for D; inferred from D;/D ) = 0.67 (E) or 0.70 (S) (de Vaucouleurs 1956b). 
(4) the D scale for E is non-linear: log D(Ia) = 0.75 log D(XIb) + 0.10. 
(5) the R scale for E is irregular: R(la) = R(XII) for R > 0.7, 
= R(XII) + 0.10 for R < 0.7. 
(6) the R scale for S is irregular: R(la) = R(XIIIb) for R > 0.6, 
; = R(XIIIb) + 0.08 for R < 0.6. 
0 A similar effect, perhaps with a different ampli- 
A tude, applies to the threshold brightness m, that 
mM, defines the visible edge of a nebular image; when 
the linear size (in millimeters on the plate) is 
large, and for a given set of observational circum- 
1 stances, m, depends only or primarily on the 
density gradient, i.e., on the luminosity law, in 
the nebular image; in other words, for a given 
galaxy type, m: is a constant, characteristic of 
2 the micrometric dimensions system which then 
0:5 1 5 10 approximates a photometric system. 


Figure 2. Variation of threshold surface brightness as a 
function of image diameter for ellipticals and spirals from 
equation (7) and for uniform images after Hubble (1932). 


On the other hand, when the linear size is small 
the threshold brightness m1 varies with image 
size as in Figure 2 so that micrometric dimensions 
measured on small scale photographs do not 
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approximate any given system of photometric 
dimensions; the measurements of small galaxies 
come out too small compared with the system 
of brightness dimensions defined by the larger 
galaxies. 

It can be shown that this effect is measured 
by the additive constant B in equation (4). For 
example, the relation between the micrometric 
dimensions D in the Heidelberg survey (series Ia, 
scale 1‘7/mm) and the brightness dimensions at 
m, = 25.0 mag/sec? are of this form: 


Dy(m,=25) =2.4D+1'2=2.4(D+0'5)  (E), 
Dy(m=25) =1.25 D+0'6=1.25 (D+0'5) (S). 


Thus the brightness diameter D, relative to a 
constant threshold m, is of the form 


Dy = (D+ 8), (5) 


where 6 is a small constant. The relative error 
in D is then given by the expression 


Alog D = log Di/D = log (1 + 6/D). (6) 


The corresponding variation of the threshold 
brightness is 


Am, = K log (1 + 6/D) (7) 


if K is the slope of the relation m, = f(log D). 
This slope can be computed for the luminosity 
laws (1) and (2) if it is assumed, from the data 
in Table I, that a ~ 1 (EZ) or 2 (So, S) for A = 1 
(see section 6). The mean slope near A = I is 
then K = 5 (£) or 9 (So, S); the values of Am, 
so computed are shown in Figure 2 where the 
curve for spirals accidentally coincides with 
Hubble’s empirical relation. The similarity of the 
relations Am, = f(D) with Hubble’s confirms the 
interpretation of the constant B in (4). 

This result makes it possible to adjust the 
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coefficients of the empirical relations between two 
systems of micrometric dimensions. If the gen- 
eral relation between brightness and micrometric 
diameters is 


Dim) =aD+b6=a(D+5) (8) 


and if the relative error in D is, for a given galaxy 
type, the same function of the linear size of the 
image on the plate for all series, 6 must be pro- 
portional to the plate scale S in minutes of arc 
per millimeter, i.e., 


5 = Sbo. (9) 


It follows that the relation between two systems 
of micrometric diameters is 


a,(D, —- S160) = a2(D» SF S250) 
or 


a a 
Dh = 7, + do (2 So — S:) (10) 


which is of the form (4). The values 6 = 0/5 for 
series Ia (S = 1/7/mm) suggest that 69 = 0/3 
for all galaxy types and this is confirmed by the 
satisfactory fit of the empirical correlograms 
given by the adjusted coefficients based on this 
value of 60, as listed in Table III. It is clear, 
however, that very small values of D, say less 
than about 36 or 1 mm on the plates, are of lower 
than average accuracy. 

The objects in Table I whose luminosity pro- 
files are known on an absolute scale have been 
used for a provisional calibration of the scale 
coefficient A. A final calibration will require more 
standards, especially for spirals, than are pres- 
ently available. The mean surface brightness 
thresholds corresponding to the micrometric ma- 
jor diameters in several series are listed in Table 
IV. Two values are given: m; is the uncorrected 


TABLE IV. MEAN THRESHOLD BRIGHTNESS OF MICROMETRIC DIAMETER SYSTEMS 


Ellipticals 
mi 


Series mi n m1 
Ia 2085 23.2 9 Danii 
Ib 22105) -23':0 4 23.90 
Ic Zeal 22.7 7 22.86 

IIb — — 25.62 

IV 22n02;) 2310 4 23.04 

VI 22528, 22.6 5 23.44 

VII Do — — 24.70 

VIII D; Z2r.1es- 9 21.3 2 22.01 
=A 24.08: 24.1 2 24.67 
IX D; 225405 22.45 6 — 
—- Do 24725. 24.25 6 


Lenticulars 
mi’ 


Mean diameter 
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Figure 3. Relation between mean threshold brightness m, (mag/sec?) and scale factor A. Empirical values (left), 
adjusted values (right). A = 1 for system la. 


average for the mean diameter listed in the last 
part of the table; m:’ = m, + Amy, is the value 
corrected for the variation Am, of threshold 
brightness as indicated by Figure 2, Both values 
appear in Figure 3 where mz; is plotted against 
the empirical value of A (left) and my’ against 
the adjusted value of A (right) as listed in Table 
Ill. The crosses are additional data from three 
series of brightness and microphotometric di- 
ameters. 

The theoretical relations m,(A), equivalent to 
log B = f(r), corresponding to luminosity laws 
(1) and (2) are shown by the curves in the right 
hand part of Figure 3. The average threshold 
brightness corresponding to A = 1 is m, = 23.0 
mag/sec? (E) or 24.2 mag/sec? (So); the latter 
may also apply to spirals, but a direct verifica- 
tion will be needed. Note that the threshold 
brightness is not the same for galaxies of different 


luminosity laws even when measured on the same, 


plates. 
Finally, if a general luminosity law charac- 


terizes ellipticals, and another characterizes spi- 
rals in their outer regions, then a statistical 
relation should exist between the scale coefficients 
A relating to each type. A plot of A(E) against 
A(S) shows that this is indeed the case and the 
mean empirical relation A(S) ~ VA(E) is in 
good agreement, at least for A 2 1, with the 
theoretical relation derived from the luminosity 
laws (1) and (2) if one assumes as before that 
a = 1 (£) or 2 (S) for A = 1. This relation was 
also used as a guide to homogenize the values of 
A in Table III. 

4. Provisional reduction to a standard system. 
A provisional reduction of series I a, b, c; Il a, b; 
Ill a, b,c, de 1V: Vis V1 a; by cee xe ie 
and XIII a, b, c to the corrected Heidelberg sys- 
tem (Ia) was made by means of the non-adjusted 
coefficients A, B of Table III. The corrected 
Heidelberg system is defined as follows: 

(i) the axial ratio R(la) is corrected according 
to Holmberg’s precepts (1946) which give the 
following relations between the photometric ratio 
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TABLE V. REDUCTION OF MICROMETRIC DIMENSIONS TO STANDARD SYSTEM 


NGC 4486, Eo-1p, mr = 11.0 


Observed Corrected 

Series D R log D —log R 
la Buea L200 0.56 0.00 

Ic Dae OOF Oy ORMIY/ 
Ila eels | SRaCeto) 0.5I 0.00 
IIb FT O.90 0.63. 0.05 
TE Be5e 51200 0.59 0.00 
x BuO 1-00: 0.47. 0.00 
XII Zen 7 LOO 0.50 0.00 
Mean 0.54 0.03 
Aver. dev. _ +05 -05 


Corrected D = 3/5 + 022 (p.e.) 
— d/D =0.93 + 0.03 (p.e.) 


b/a and the micrometric ratio B/a as measured 
by Reinmuth (1926) 


bla Bla (E) Bla (S) b/a Bla (E) Bla (S) 
1.0 1.0 1.0 0.5 0.397 0.429 
0.9 0.873 0.882 0.4 0.289 0.322 
0.8 0.750 0.766 Ona 0.192 0.223 
0.7 0.630 0.652 On2 OphUg On137 
0.6 0.512 0.540 O.1 0.049 0.062 


A comparison of the values so corrected with 
photometric determinations for the objects of 
Table I confirms the validity of Holmberg’s 
corrections ; 

(ii) the diameter D(Ia) is corrected by +0/5 
according to the precepts of section 3 to give 
equivalent brightness dimensions at m1 23.0 
mag/sec” (E) or 24.2 (S). In order to reduce to 
a common threshold brightness, say m,; = 23.5 
mag/sec”, corrected diameters should be multi- 
plied by about 1.25 (E) or 0.80 (S); reduction 
factors to m,; = 25.0 mag/sec? would be 2.4 (E) 
or 1.25 (S), but the application of such large 
corrections is not advisable. 

Mean corrected dimensions were computed for 
about 1200 Shapley-Ames galaxies listed in two 
or more series. Table V gives two examples, for 
an elliptical and a spiral. The complete list of 
corrected diameters for over 2000 galaxies based 
on all available series will be published later. 

Probable errors computed for all objects whose 
means are based on at least four values indicate 
that the individual accuracy of dimensions data 
varies relatively little from one series to another ; 
it is generally within +30 per cent of the mean 
p.e. = 0.047 in log D/d, or +12 per cent in d/D, 
independent of type; it is within +50 per cent 
of the mean p.e. = 0.061 in log D, or 15 per cent 
in D, for ellipticals, and of the mean p.e. = 0.040 
in log D, or 10 per cent in D, for spirals. The 
different accuracy of diameter data for spirals 
and ellipticals is a consequence of the different 
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NGC 4216, SAB(r:)b, mr = 10.8 
Observed Corrected 

Series D R log D —log R 
la ae Ory 0.90 0.70 

Ib Fa. Ogi) Ox908) 0-74 

Ic 6'5 0.08 0.88 0.86 
IIb 9/0 0.20 0.83. 0.70 
IV Seen meOmAs 0.96: 0.69: 
x Omen O SL 0.89 0.68 
Xla COLON OnT7: 0.81 0.64 
XII Hag) Onitss 0.90 0.67 
XIle ioe Oe20 0.91 0.57 
Mean 0.89 0.69 
Aver. dev. _ -03 .05 


Corrected D = 7/7 + 0/2 (p.e.) 
—— 1 d/ DF 020-0701 (pies) 


luminosity laws; the corresponding common value 
of the p.e. of the threshold brightness m, is about 
0.4 mag. 

Thus, if independent dimension data are avail- 
able for a bright galaxy from, say, four different 
series, the mean corrected diameter is known 
with a p.e. of 5 or 7 per cent (S or E) and the 
axial ratio with a p.e. of 6 per cent (S, E). Fur- 
ther, if adequate calibration data are available, 
the surface brightness of the corresponding iso- 
phote is known with a p.e. of 0.2 mag. This 
result encourages the belief that micrometric 
dimension data so corrected could be used as 
distance indicators (see section 5). 

In preparation for this application some factors 
affecting the apparent diameters of galaxies must 
be considered. Figure 4 shows the correlation 
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Figure 4. Relation between mean apparent diameter 
and axial ratio for galaxies of total magnitudes 11.0 < mr 
eo. T OO: 
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between mean apparent diameter D in the cor- 
rected Heidelberg system and axial ratio for 
objects of types E, So-Sa, Sb, c, d and total 
magnitudes 11.0 < mr < 12.0; the relation is 
essentially the same for all types, 


log D = log D(o) + 0.40 log D/d_ = (11) 


D(o) is the face-on diameter, i.e. the diameter 
that would be observed for D/d = 1 orlog D/d = o. 
This relation agrees closely with that derived 
from fewer data in the Mount Stromlo survey 
(de Vaucouleurs 1956b) ; on the average a galaxy 
seen edgewise with an axial ratio R = d/D = 0.2 
appears with a maximum diameter almost twice 
its face-on diameter D(o). This, of course, is the 
cause of the apparent excess of edgewise objects 
in many statistical surveys (Wyatt & Brown 
1955). 

Figure 5 shows the correlation between appar- 
ent diameter and galactic latitude for objects of 
axial ratio R > 0.5 and total magnitudes 10.5 < 
mr < 12.0. A variation of the order of A log D = 
+0.15, or 40 per cent in D, between 6 = 20° and 
b = 70° is indicated, and confirms qualitatively 
earlier results of Reiz (1941) based on the un- 
corrected Shapley-Ames magnitudes and diam- 
eters. The mechanism of this variation has been 
discussed earlier (de Vaucouleurs 1957a), but the 
exact form of the relation remains to be de- 
termined. 


60° 


0 30° 90° 


Figure 5. Relation between mean apparent diameter and 
galactic latitude for galaxies of axial ratio R > 0.5 and 
total magnitude 10.5 < mz < 12.0. 


5. Apparent diameters as distance indicators. 
In high galactic latitudes the relation between 
apparent diameter and total magnitude can be 
written 


mr +'5 log D(o) = const., (12) 
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where the constant depends only on type and 
surface brightness. 

The scattered distribution of the S—A galaxies: 
as a function of corrected Harvard magnitude 
m, and of conventional surface brightness m,’ = 
m, + 2.5 log Dd was illustrated previously (de 
Vaucouleurs 1956a, 1957a) for the dimension 
data of the S-A catalogue. A conspicuous reduc- 
tion of the scatter takes place when the mean 
corrected diameters are used as shown by Fig- 
ure 6 giving plots of mr vs. D(o) for galaxies in 
the 12 degrees core of the Virgo cluster (6 = 
+75°). The constant of relation (12) is 13.0 for 
ellipticals, 13.5 for So-Sa and 13.8 for Sb, c, d; 
for the last group the dispersion of individual. 
values hardly exceeds 0.2 mag. 

This result does not imply a near uniformity 
of intrinsic surface brightness among galaxies of 
a given type; it arises simply from the fact that 
brightness dimensions are directly correlated and 


0 05 1 


Figure 6. Correlation between total magnitudes and 
“face-on”’ diameters of galaxies in the Virgo Cluster. 
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total magnitudes are inversely correlated with 
surface brightness, so that the opposite variations 
of the two quantities balance out to a large 
extent. It is of some interest that for the bright- 
ness dimensions at m; = 25.0 mag/sec? the con- 
stant would be about 14.0 for all galaxy types 
(excluding low luminosity elliptical and irregular 
dwarfs); further, since the mean absolute mag- 
nitude of galaxies of all types is also a constant 
(de Vaucouleurs 1958), it follows that the mean 
brightness dimensions at m; = 25.0 mag/sec? are 
also about the same for all types (with the same 
exclusions). 

Because of the close correlation between D(0) 
and mr either can be used as a distance indicator. 
Here the speed and simplicity of determination 
of micrometric dimensions is an important ad- 
vantage in comparison with the slow and labori- 
ous photometry. 

The practical possibility of using reduced diam- 
eters as statistical distance indicators was first 
explored in the Mount Stromlo Survey (de Vau- 
couleurs 1956b, pp. 29-32) ; the relative distances 
of seven clouds or regions in the southern sky so 
determined were placed on an arbitrary distance 
scale. The provisional mean diameters now avail- 
able for nearly all S—A galaxies in both hemi- 
spheres provides for a tie with the current distance 
scale (Sandage 1958) which is essentially based 
on northern hemisphere objects. Table VI gives 
the mean reduced diameters and total magni- 
tudes of the largest spirals in some nearby groups 
whose distance moduli are known; no distinction 
between apparent and corrected moduli needs to 
be made at the present degree of accuracy, and 
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with this accuracy the rounded-off mean values 
M (pg) = —20.0 and D(o) = 20 kpc may be 
adopted. Here again D is the mean diameter in 
the corrected Heidelberg system and refers to 
m, ™ 24.2 mag/sec?. Table Vla gives conversely 
the distance moduli derived with these values 
for M31-M33, the Ursa Major group, the Virgo 
Cluster, and for the compact cloud of spirals in 
Grus (region VII in the Mount Stromlo Survey). 
The mean modulus m — M = 31.2 places this 
cloud at a distance of 17.5 Mpc on the current 
scale, as compared with 7.0 Mpc on the distance 
scale used in 1956. It follows that the mean dis- 
tance of region VI (Pavo-Indus) which was used 
as a unit equal to 10 Mpc in the Mount Stromlo 
Survey is now 26.5 Mpc and the adopted dis- 
tances of the other groups are similarly changed ; 
in particular the estimated distance of the Fornax 
cluster is now about 8 Mpc. More accurate values 
of the relative distances of these and other group- 
ings of bright galaxies will be derived from the 
revised dimension data in the final catalogue. 

6. Asymptotic magnitudes. The initial impetus 
to secure better diameter data was the need to 
improve the extrapolation correction m, — mr 
that must be applied to integrated magnitudes 
m, to determine the total (or asymptotic) mag- 
nitudes mr. 

If the luminosity laws (1) and (2) apply the 
correction is readily computed. Consider for sim- 
plicity circular objects. For ellipticals the total 
luminosity is, by integration of (1), 


Lr = on [ BG@)rdr = 7.215 7Beré (E) (3) 
0 


TABLE VI. MEAN DIAMETERS AND MAGNITUDES OF LARGE SPIRALS 


Group D(o) mr mr + 5 log D(o) 
Local 88/5: 5.0: TAGS: 
Ursa Maj. m3 8.4 14.4 
Virgo 5/0 10.4 13.9 
Grus BES ee) 1307, 

Remarks: 


(1) mean of M31, M33. 
(2) mean of M51, M81, Mror. 


m—M M D (kpc) Remarks 
24.5 —I19.5: 20.5: (1) 
28.0 —19.6 20 (2) 
30.5 —20.1 18 (3) 
se = = (4) 


(3) mean of Io largest (4’ to 6‘5) or brightest (9.7 to 11.0) out of 36 in core of cluster. 
(4) mean of 5 largest (3/0 to 5/5) or brightest (10.7 to 11.3) out of I9 in region VII of Mt. Stromlo survey. 


TABLE VIa. DISTANCE MODULI FOR M = —20.0 AND D = 20 KPC 
Group Local UMa Virgo Grus 
From mr 25.0: 28.4 30.4 31.0 
From D(o) 24.4: 28.0 30.7 31.4 
Mean m — M 247s 28.2 30.6 21.2 
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and the fraction L,/Ly emitted within the iso- 
phote of radius 7 is 


Bi) Ley Le 


RG ace Bas 
= I -e+(1 a iene ee +=) (14) 
if x = 7.668 a!/4; by definition r = 7. for k = 3, 
ANGUBei— 9D (fe) f 
For spirals the total luminosity of the expo-_ 


nential component is, by integration of (2), ; 

(S) Lr = 3.803 mBeré (15)f 

0 and the fraction emitted within r is | 
-05 0-0 05 log x 1.0 

Rr) =1—-U+n)e7 (16) 


Figure 7. Relative integrated luminosity me — mr = : 4 
—2.5 log k(a) as a function of reduced diameter a = D/D, if x = 1.6785 a and, as above, k = 3 fora = r/r, 
for the luminosity laws (1) and (2). D, = effective di- — 1, Figure 7 gives the variations of R(r) and of 
ameter. 


mM, — mr = —2.5 log kf) (17) 


(a) 
-3 
-05 0 logA-logD 05 -0:3 0 log A/D 05 


Figure 8. (a) Mean slope Am,/A log A of integrated luminosity curve of Sb, Sc spirals as a function of diaphragm 
aperture A = 2r in the data of Pettit (1954). (b) Empirical incompleteness correction curve m, — mr as a function 


of log A/D. D is the mean micrometric diameter in the corrected Heidelberg system (D & 2.6 D.). 
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TABLE VII. NGC 157, Sc, D = 3/0 


A = 29 Bas All 
My 11.19 Il.O1 
log A/D —0.115 +0.135 
mM, — mr 0.17 0.03 
Moorr II.02 10.98 


as a function of @ for the two cases. The first, 


can in ellipticals, gives k(r) > 0.90, 

— mr < 0.10 mag. for a > 5, R(r) > 0.97, 
a — mr < 0.03 mag, for a > 10; the second, 
dominating in late-type spirals, gives k(r) > 0.90 
for a > 2.5, k(r) > 0.97 for a > 3. Because of 
_ the rapid convergence of the luminosity laws (1) 
and (2) the extrapolation ‘‘to infinity”’ involved 
in the definition of the asymptotic magnitudes 


can in fact be limited to more realistic distances. 


In order to apply Figure 7 .to the statistical 
Borrection of measured integrated magnitudes 
the mean ratio between effective and micro- 
metric dimensions (in the corrected Heidelberg 
system) is needed. A direct comparison for ob- 
jects in Table I gives D ~ 1.3 D, for ellipticals, 


and D~ 1.7 D, for lenticular galaxies. For spi- 


rals an indirect determination was made by 
means of the photoelectric magnitudes m, meas- 
ured by Pettit (1954); Figure 8a shows the 
variation of Am,/A log A asa function of log A/D 


_ for Sb, Sc spirals measured through at least two 


different apertures A = 27; only objects of axial 
ratio R = d/D > 0.5 and galactic latitude |b| > 
45° are included. 

By integration the empirical relation m, — mr 
= f(log A/D) shown in Figure 8, 0 is obtained; 
it leads to D ~ 2.6 D, and consequently log a = 
log A/D + 0.42; with this value the theoretical 


relation (S) of Figure 7 coincides almost exactly 


with the empirical relation of Figure 8, b. This 
result confirms the validity of equation (2), at 
least in the intermediate and outer regions of 
spirals. 

Table VII gives an example of determination 
of the asymptotic magnitude of a spiral meas- 
ured by Pettit (1954). Asymptotic magnitudes 
so determined will be given in the reference cata- 
logue of bright galaxies now in preparation. 

I am much indebted to Mrs. A. de Vaucouleurs 


BT) Adopted 
10.98 
+0.28 

0.01 

10.97 mrp = 10.99 + 0.02 


for assistance in the numerical work and to the 
many persons who over the past ten years have 
made possible the collection of basic data at 
several observatories. Details of these data will 
be reported separately. 
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Abstract. A multiple-exposure technique, employing separated photographic and photovisual exposures on the same 


plate, has been tested for use in galactic cluster photometry. Results, to a given accuracy in B— 


V, are achieved with © 


less work than that required for conventional single exposures. 
A special emulsion is recommended, and photometric properties of the Naval Observatory’s 40-inch reflector are” 


discussed. 


1. Introduction. A photometric investigation 


of a large number of galactic clusters is being 


carried on jointly by the Lowell and U. S. Naval 
Observatories. In the course of this program, we 
intend to measure magnitudes and colors of some 
20,000 stars on plates taken with the Naval 
Observatory’s 40-inch reflector. Because of the 
size of this program, we have sought a method of 
photographic photometry which will be as eff- 
cient as possible. This report deals with a method 
of multiple exposures on a single plate. Haro has 
used this method qualitatively with great effec- 
tiveness (Haro and Herbig 1955, Haro 1956). 
Quantitatively, this technique of measuring colors 
has two advantages over the conventional method 
of single exposures; first, the number of identi- 
fication settings in photometering the plates is 
reduced by a factor of two, and secondly, differ- 
ential emulsion and processing errors are reduced. 
Our multiple exposure plates are obtained in the 
following way: Plates of a given field are taken 
in pairs. The first exposure is made through a 
Wratten 47 plus Schott GG13 filter combination 
on a D sensitized emulsion. The plate is shifted 
0.5 mm (15’”), and a second exposure is made 
through a Schott GGr1 filter. The second plate 
of the pair is then taken with the, exposures in 
reverse order. The relative exposure times are 
adjusted to give approximately equivalent images 
in blue and yellow for the maximum of the fre- 
quency distribution in color for the field.. Images 
are evaluated with an Eichner photometer (Eich- 
ner ef al. 1947; Cameron 1951), and by means of 
distributed photoelectric standards. 

2. Selection of an emulsion. In order to appraise 
the multiple-exposure method, we obtained a 
series of test plates with various emulsions. Each 
test plate was multiply exposed to a selected 
field with two equal exposures separated by 0.5 
mm (15’’). Pairs of images of a magnitude 
sequence selected within 25 mm of the plate 
center were then measured with the Ejichner 
photometer. Each image of each pair was meas- 


{ 


if 


ured twice so that_our measuring error could be, 
evaluated. The calibration curves for the various — 


emulsions were all quite similar, giving a scale 


of about 2 magnitudes per 1000 divisions of the 


iris scale. 

Results of these double-exposure tests can be 
presented by plotting the mean Eichner reading, 
R, against AR, the difference in Eichner readings, 


for the various pairs in a sequence. For a perfect. 


- plate with identical exposures, AR would be zero 


in every case. In practice, the R vs. AR diagrams 
define a line about which there is some scatter. 
The scatter of the points may be used as an index 
of precision, the slope of the line depending on 


the relative circumstances of the two exposures — 


such as seeing, focus, and guiding. The results 


of several double exposure tests are summarized | 


in Table I. 
TABLE I. RESULTS OF DOUBLE EXPOSURE TESTS 
Exposure 
No. of p.e. one p.e. one for given 
Emulsion plates measure difference limit 
103a—O 2 +o™01I £07033 I 
Ila—O 2 0.010 0.021 2 
III-G I 0.011 0.021 
IV-D I 0.012 0.04 
Separation 
Neg. Type I D 0.015 0.032 
33 6 0.010 0.016 6 
Exp. D 5 0.010 0.016 6 


Sample tests of the kind leading to the results 
shown in Table I are illustrated in Figures I 
and 2. Our initial tests showed that the Eastman 
33 plate was superior to normally available spec- 
troscopic emulsions for this work. Accordingly, 
we asked Mr. William Swann of the Eastman 
Kodak Company if D-sensitized 33 plates could 
be supplied for use in multiple exposure photog- 
raphy with blue and yellow filters. This emulsion, 
designated Experimental D in Table I and Fig- 
ure 2, was supplied by the courtesy of the East- 
man Kodak Company. Additional material of 
this type, which has been obtained for the cluster 
program, is designated Experimental CC—D. 


| 
q 

i! 
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Figure 1. Eichner measures for sample double exposure 
test cf 103a—O emulsion. 


Further experiments have shown that the photo- 
metric properties of this emulsion are similar to 
those of the 33 plate. 

The last column of Table I gives the relative 
speeds of four emulsions of interest. It appears 
that the Experimental emulsion is inconveniently 
slow. This emulsion, however, responds well to 
sensitization by baking. We have realized gains 
of one magnitude for exposures as short as I5 


minutes when these plates have been baked for 
24 hours at 65° C. This gain in speed is achieved 


without deleterious effects on fog level or grain 
size. Response to baking, and consideration of 
the fact that the ratio of exposure time to meas- 
uring time is small, are factors which make the 
use of the Experimental emulsion practicable. 

We found that the index of precision of the 
Experimental emulsion is not improved by the 
use of fine grain developer as compared to devel- 
opment in D-19. 

Double exposure tests made with either the 
Wr. 47 + GG13 or GG11 filters used in the pro- 
gram work show an increase in the probable error 
of the difference between images of a pair, the 
error introduced by the filters being something 
of the order of +o0™o10. 

3. Photometric properties of the Ritchey-Chrétien 
reflector. The 40-inch reflector has been described 


Mme ASTRONOMICAL JOURNAL 


411 
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background 


AR 


Figure 2. Eichner measures for sample double exposure 
test of Experimental D emulsion. 


by Hall and Mikesell (1950), and prior to that 
by Ritchey (1928) and others. This Cassegrain- 
like telescope has a focal surface concave towards 
the secondary with a radius of 65 inches, and a 
scale of approximately 30’’/mm. By bending thin 
plates in a vacuum back, fairly good images can 
be obtained over a 1°5 field. However, we prefer 
to use flat plates for convenience. 

When flat plates are used, the field error attrib- 
utable to curvature of field is large in comparison 
to the geometrical vignetting of the optical sys- 
tem. The photometric errors produced by this 
curvature of field were determined by Eichner 
measures of well distributed photoelectric stand- 
ards on focus plates. The focus plates were taken 
in such a way as to provide 10 images covering 
the range of focus present within a 50 mm (25’) 
radius from the center of a flat plate. The error 
to be expected from curvature of field varies 
according to the relative brightness of the star 
measured. Figure 3 shows errors determined for 
images between I and 4 magnitudes above the 
plate limit as a function of the distance from the 
center of the plate. These values obtain when 
the best focus is established by knife edge meas- 
ure at the center of the plate. Note that unsatu- 
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unsaturated images 


saturated images 


(s) 5' 10° 15' 20' 


te) 10 20 30 


Figure 3. Photometric errors due to curvature of field 
as a function of distance from the center of the plate. 
Crosses, open circles, and closed circles represent mean 


40 mm. 


data for images I, 2, and 4 magnitudes above the plate . 


limit, respectively. 


rated images measure fainter at large distances 
from the plate center when defocused, and satu- 
rated images measure brighter. In order to avoid 
scale errors, the measured field must be kept 
within a diameter of 20-25’. This field is ade- 
quate for the present photometric program. 

Astronomical seeing, insofar as aperture, focal 

length, and site are concerned, may be regarded 
as a photometric property of the telescope. 
“‘Seeing’”’ is measured with the same knife edge 
apparatus that is used for focusing by the Fou- 
cault test. The position of the knife edge is regis- 
tered by a dial indicator. A K.E. cut, i.e., the 
distance travelled by the knife edge from the 
position where the first shadows appear on the 
mirror to the position where the last bright spots 
disappear, is an index of the performance of the 
telescope as limited by the seeing and the figure 
of the mirror. This Foucault shadow measure 
gives more consistent results than direct meas- 
ures of the diameters of images. Furthermore, 
_K.E. cuts measured in this way give uniform 
results over large ranges of magnitude or tele- 
scope aperture provided that the K.E. cut is 
made slowly enough to integrate all wave front 
disturbances at the mirror. An effective focal 
length of 5 meters or more is required for con- 
venient application of this method. 

On a given night, the K.E. cut will vary with 
the air mass at which the observation is made. 
Examples of this variation are shown in Figure 4 
for nights of ‘“‘good”’ and ‘“‘bad”’ seeing. In order 
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Figure 4. Variation in seeing with change in air mass for ; 
two nights with different seeing conditions. 


to provide a uniform basis for night-to-night com-_ 
parisons, K.E. cut measures are reduced to the 
zenith by the following relation: 


(K.E.)z = (K.E.)a.m. — (A.M. — 1)0.05. 


The K.E. cuts are expressed in millimeters and _ 
A.M. refers to the air mass at which the observa- 
tion was made. 

A frequency distribution of K.E. measures for 
100 nights distributed more or less at random 
over a three year period is shown in Figure 5. 
The three nights illustrated at the tail of the 
distribution show how miserable conditions can 


100 nights 


[oH] 0.2 0.3 o4 0.5 0.6 


K.E. cut mm. 


Figure 5. Frequency distribution of seeing measures 
for a three year period. 
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Figure 6. Percentage of starlight admitted through various 
sized apertures on two different nights. 


get after storm front passages. In order to relate 
these data to terms more meaningful than units 
of K.E. cut, photoelectric calibrations of these 
seeing measures have been made in two ways. 
One approach has been to measure a bright star 
with a series of different sized apertures. Mean 
results of series of such measures for two nights 
are shown in Figure 6. These data show that an 
aperture diameter equal to the K.E. cut on a 
particular night will admit approximately 90% 
of a star’s light. 

Photometer tracings of actual knife edge cuts 
have provided another means of calibration. 
From these tracings, the knife edge excursion 
required to occult the central 68% of the light 
of an image can be determined. These measures 
are equivalent to those reported by Meinel (1958) 
in which occulting bars were used to assess the 
sizes of images. If the same correction for instru- 
mental scattering applied by Meinel is made, it 
is found that: 


oe “Se SPOKE.) sam 


o”’ is the width of an occulting bar in seconds of 
are which will obstruct 68% of the light of an 
image if correction is made for instrumental 
scattering. 

If we apply this latter calibration conversion 
to the data illustrated in Figure 5, we find that 
the average o”’ for our site at Flagstaff is 1/4. 
This same value has been reported by Meinel 
for 150 nights distributed over one year’s time 
at Kitt Peak for observations of Polaris. 

The actual photographic performance of the 
telescope as related to seeing is illustrated in 
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Figure 7. Diameters of photographic images a few tenths 
of a magnitude above the plate limit as a function of seeing 
for 100 103a—O chart plates. 


Figure 7. Diameters of photographic images a 
few tenths of a magnitude above the plate limit 
are plotted against the K.. cut index of seeing 
for 100 103a-O chart plates. Diameter measures 
were made with a 60-power binocular microscope 
so as to include all grains attributable to the 
images. Each point on the diagram represents 
the mean of several such measures for one plate. 
The large scatter in the diagram can be explained 
by errors in focus and guiding, and variations in 
seeing during a given exposure. Because there is 
not a one-to-one correspondence between knife- 
edge cut measures and the limiting image sizes, 
it is apparent that only the intensity maximum 
for faint images is effective in producing a photo- 
graphic image. The “‘washing-out”’ of the image 
intensity profiles as a function of seeing quality 
results in loss of limiting magnitude for a given 
exposure because of the relative reduction of 
central intensity and effective increase in area 
of the integrated optical images. Knife-edge cut 
measures are of some use in anticipating the ex- 
posure required for a given limiting magnitude. 

4. Preliminary results of multiple exposure pho- 
tometry. Preliminary results of photoelectrically 
calibrated multiple-exposure photometry for 20 
galactic clusters are presently available for evalu- 
ation. External errors for the photographic work 
have been estimated by comparing photographic 
and photoelectric results for about 30 stars in 
each cluster. Comparisons are made as follows: 
Separate calibration curves are constructed for 
the means of the photographic and photovisual 
measures for pairs of plates. Instrumental mag- 
nitudes, 6 and v, may then be read from these 
curves. We have so far found that: 


Sv, 


and 
B-—-V=4,+A42(b —»), 
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where A, is a function of the photovisual Eichner 
photometer reading, and 


A» = 0.95(I + 0.03 A.M.). 


The term in the parentheses in the latter equa- 
tion takes into account differential extinction as 


a function of the mean air mass, A.M., and the 
scale factor, 0.95, has been determined from re- 
sults for 300 stars. Further details of the reduc- 
tion procedures will be given when the final 
results are published. 

Probable errors are given in Table II for results 
from pairs of plates. Arp (1955) has pointed out 
that Ila-O plates give better results than 103a-O 
plates. Our data show that a further improve- 
ment can be realized by the use of the Experi- 
mental emulsion. These data also show that the 
multiple exposure technique gives better results 
than an equivalent amount of material on sepa- 
rated single exposures since the errors in B—V 
are somewhat less than those in V rather than 
being V2 times as great. We also conclude that, 
for a given amount of plate material, photometric 
errors cannot be substantially further reduced by 
improvement in emulsions as the Experimental 
emulsion error is small in comparison to the total 
error introduced by other steps in the photo- 
metric process. 


TABLE Il. EXTERNAL PROBABLE ERRORS FOR MULTIPLE 
EXPOSURE PLATE PAIRS 


No 
Emulsion pairs p.e. (V) p.e. (B—V) 
103a—D 2 =-0™036 07033 
Ila—D 16 0.031 0.026 
Exp. CC-D 2 0.027 0.022 


5. Summary. We have found that a multiple- 
exposure method of measuring colors photo- 
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graphically gives better results than an woudl 
number of separated exposures. This is presum- i 
ably brought about by a reduction in the effects 
of emulsion and processing irregularities. The 
method is also more efficient in that identifica- 
tion settings with the plate photometer are cut 
in half, for B— V measures, as compared to meas- 
urement of single exposure plates. 

Weare now using the multiple-exposure methoml 
for three-color observations in our galactic cluster 
program. ; 
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and secondary components respectively. 


Introduction. More than one revolution has 
been completed since the binary character of 
‘10 Ursae Majoris, 854™2, +42°11’ (1900) was 
‘discovered by Kuiper (1935). The visual magni- 
‘tudes of the components given by Eggen (1956) 
are 4.19 and 6.02; the spectrum has been classi- 
fied as F5V by Johnson and Morgan (1953). The 
recently determined orbit by Baize (1955), P = 
22%20,a = "61, combined with the present Sproul 
‘material permits a determination of the masses 
of the components with substantially greater 
‘accuracy than was hitherto possible. The present 
‘determination includes the Sproul material used 
in the earlier solutions for parallax (Miller 1918) 
and for parallax and mass ratio (van de Kamp 
1947). 
_ Material and measures. The procedure described 
‘in previous Sproul papers has been followed (van 
‘de Kamp and Lippincott 1949). The material 
includes 608 exposures on 186 plates, taken on 
67 nights, resulting in a total weight of 142. The 
‘material is summarized in Table I. The regular 
exposure time was 10 minutes for the interval 
1912-1916, 2 minutes for the interval 1939-1945 
and 1% minutes since 1945. The dependences 
and their yearly changes are given in Table II 
which also includes the standard frame and posi- 
tion of 10 UMa for 1940.0. The spectra of the 
reference stars were kindly furnished by Dr. 
Vyssotsky. The brightness of 10 UMa was re- 
duced to match the average magnitude 10.8 of 
the reference stars by the use of a rotating sector 
of .20 per cent opening. The plates were meas- 
ured on the St. Clair-Kasten long-screw meas- 
uring machine in the directions of right-ascension 
(x) and declination (y) for the approximate 
equator of the year 2000. The measurements 
and plate reductions were made by Mrs. John E. 
Houtman for the material from 1912-1947, and 
by Mrs. Mary Jackson for the material 1949- 
1959. Mrs. Jackson also remeasured the plates 
from the 1940-43 seasons. No significant differ- 
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PARALLAX AND MASS RATIO OF 10 URSAE MAJORIS 


| FROM PHOTOGRAPHS TAKEN WITH THE 24-INCH SPROUL REFRACTOR 


By SARAH LEE LIPPINCOTT 
| Sproul Observatory, Swarthmore College, Swarthmore, Pa. 


| Received May 1, 1959 


Abstract. Measurement and reduction of photographic plates taken with the Sproul refractor over the interval r9r2- 
|1959 yield +0721 + "0047 (p.e.) for the relative parallax and +1702 + "oo40 for the semi-axis major of the photo- 
centric orbit. The adopted value for the absolute parallax, +7074 + "003, leads to masses .63 © and .49 © for the primary 


ences between the two sets of measurements were 
found; straight means are used with no increase 
in weight. 

Solution for parallax and photocentric semi-axis 
major. The separation of the components never 
exceeds 77 or about .o4 mm; the measured posi- 
tion of the blended image has been assumed to 
represent the center of the light intensity or 
photocenter of the components. The night means 
are represented by the following equations of 
conditions : 


Deen ah | meg | a0 
Y=co,+ mt + rPs + aQ;, 


where the symbols have the usual meaning; ¢ is 
counted from 1940.00. The orbital factors were 
computed from Baize’s elements (1955): 


P=) 22920 = 134-8 
6= 0.17 (i) as Deere 
T = 1949.97 Q5= 52675. 


Separate solutions in right ascension and declina- 
tion yield the following results: 


p.e. Weight 
Cx 333333 
Cy +6.57652 
Me —0.022420 = —"4231 +0003 21972 
My —0.0I1217I = —.2207 + .0002 21060 
Tx +0.003833 = +.0723 + .0057 66.05 
Ty +0.003906 = +.0737 = 0125 Tsay 
Oy +0.008558 = +.1615 —E .0069 44.74 
Qy +0.009394 = +.1773 + .0046 55-38 
p.e.I, 20.0025 = +.046 
peIy + .0018 = +.034 


A combined solution for 7 and a gives the follow- 
ing results: 


p.e. Weight 
Cx +3 -33336 
by +6.57640 
Mz — .022416 = —"42300 +"00027 22330 
My — .012166 = —.22958 +.00027 22368 
T + .003820 = +.072I =k .0047 73.65 
a + .009022 = +.1702 =k .0040 100.89 
p.e. I + .002I1. = +.040 
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Epoch 
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46. 
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.248 
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148 
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53- 

.856 

.916 
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122 
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121 
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+274 
.872 
56. 
-O16 
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TABLE I. OBSERVING DATA, MEASURED POSITIONS AND RESIDUALS 
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.87 
.88 
. 88 
.88 


.88 
.92 
OI 


x Y 
unit .0oor mm 
+3 9399 +6 8977 

9397 8961 
9391 9006 
9382 8955 
9365 9005 
+3 9354 +6 9003 
9250 8919 © 
g106 8872 
8847 8728 
___ 8809 8651 
+3 8792 +6 8700 
8842 8656 
8726 8612 
8722 8616 
8644 8587 
+3 8631 +6 8596 
3542 5809 
3371 5740 
3412 5684 
3187 5558 
+3 3106 +6 5552 
2958 5459 
2616 5341 
2627 5337 
2560 5257 
+3 2417 +6 5236 
2320 5160 
2200 5158 
2092 5049 
1810 4971 
+3 1748 +6 4935 
1738 4975 
1727 4945 
1766 4950 
1707 4922 
+3 1656 +6 4829 
1506 4868 
1433 4759 
1154 4650 
1091 4644 
+3 0862 +6 4501 
0536 4355 
0342 4191 
0216 4107 
0193 4069 
+3 0060 +6 4054 
0079 4051 — 
2 9854 3887 
9928 3944 
9820 3903 
+2 9795 +6 3897 
9812 3887 
9759 3775 
9678 3758 
9666 3781 
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TABLE I (continued) 
H.A, No., of Xx ¥ Vz Vy 
Date min. pl. exp. Obs. Pa Pé Qa 05 unit .00o1 mm unit .0oor mm p 
Apr. I —i1 4 14 Er —.83 +.43 —.73 — .68 ‘+29598 163743 +5 +1 3 
Boros 7 4 14 Po —.86 +.42 —.73 = .68 9600 3743 +10, 33 « 3 
WNovetss—258 94 16 Po -+.93 —.17 —.60 — .80 9494 3619 7/2 LA 
Pebr17 —30 4 14 Jo —.28 +.47 =—.66 — .84 9391 3628 33 ai tOmmnS 
Mewes —10 2 8 Br =—.46 -+.49° —.66 — .85 9447 3593 —-39 —24 I 
Mar. 7 +8 4 14 Br —.54 +.49 —.54 — .99 +29220 +63465 +31 —16 2 
Mate2290.—28)52)-5 Sc —.73 +.47 =—.53 —1.00 9183 3491 ae @) spits 
Dect Set 4 Jo. -+.73 -+-.07 =.42 —=1-06 9070 3372 =O ae 7 tl 
Feb. 20 —30 4 16 Wd —.32 +.48 —.39 —1.08 8978 3331 —18 =25 3 
Feb. 21 —35 4 16 Jo —.33 +.48 —.39 —1.08 8983 3342 —12 -I4 3 
Feb. 22 —35 4 16 Jo —.35 +.48 —.39 —1.08 +2 8988 +6 3351 —-6—5 3 
Feb. 26 0 4 8 Wy —.4r +.48 —.39 —1.08 8948 3345 —40 —9 2 
a = Samuel G. Barton Ho = Louis N. Howard = Joseph D. Rutledge 
= Leendert Binnendijk Jo =F. Jerrold Josties Ss = T. Paul Schultz 
i = Edwin V. Bishop Ma = Walter A. Matos S = Hannah B. Steele 
Br = Sheila V. Brown M = John A. Miller s = Walter A. Strauss 
Da = John J. Davis, Jr. N = Edward P. Neuburg T = Armstrong Thomas 
D = Roy W. Delaplaine P = John H. Pitman We H. John Wood 
fl = Edith Flather Po = William Poole, Jr. = Arne A. Wyller 
Fr = Laurence W. Fredrick Ra = F. Duane Ray “y = George B. Yntema 
Ro = Karl Hans Roth 
TABLE II. REFERENCE STARS 
Dep. 
Mpv Spectrum Diameter Xa Ys 1940 AD/yr 
10.8’ F5 137 —58.55 +12.16 363 + .00009 
ibedt Go LL, ap Bn 3240 203 + .00025 
II.0 Ko .125 +55.25 +20.24 -434 — .00034 
UMa (10.8) F5V pay. + 3.34 + 6.57 


Table I gives the residuals from the combined 
solution for each night. 

Discussion. Table III gives the mean epoch 
for the normal places of residuals in x and in y, 
together with the total weight and number of 
nights. These residuals in seconds of arc are 
plotted with respect to the photocentric orbit in 
the central portion of Figure 1. In the outer 
part of the figure, the visual observations are 
plotted with Baize’s orbit for comparison. Be- 
tween 1940 and 1945 and again between 1952 
and 1955, the photographic positions fall behind 
the orbital ephemeris. The visual material is 
stronger than the photographic material and 
shows no similar anomalies. Visual positions since 
Baize’s orbit determination have been added. 
Although all the observations near apastron tend 
to fall outside the orbit, it is probably premature 
to revise the orbit at this time. 

The present value for the relative parallax 
+0721 + "0047, supersedes the earlier Sproul 
determination based on part of the present inves- 


tigation. From the average magnitude 10.8 and 
spectrum G2 of the reference stars, at galactic 
latitude +43°, we find +004 for the mean 
parallax of the reference system (Vyssotsky and 
Williams 1948), (Binnendijk 1943). This leads 
to the value +076 for the Sproul absolute par- 
allax of 10 UMa. 


TABLE III. NORMAL PLACES OF RESIDUALS 


Vy 


Epoch =p =n unit .0oor mm 

1912.98 7 7 am os Co) 
15.48 9 9 sp dt Fino 
40.22 Ii 5 = 3 — 6 
43.28 13 5 — +10) a 
45-37 9 3 = 2 mele 
47.18 14 6 Sauls ="9 
48.37 9 3 +25 aaa 4 
50-32 9 3 =a) +20 
52.66 6 2 +10 + 8 
54-09 10 4 +14 +15 
55-20 12 5 +10 Ss a 
56.12 II 5 + 2 Oo 
57-38 10 5 a} — 6 
59-13 12 5 —16 S10; 
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Figure 1. Visual and photographic observations of 10 
Ursae Majoris. 

The outer part shows the visual observations (B with 
respect to A at the origin) from 1935 to 1957 grouped 
in yearly means, and the relative orbit (Baize 1955). 
The central part shows normal places of the Sproul 
photographic positions, over the interval 1912-1959 (Table 
III), after eliminating parallax and proper motion; the 
center of mass of the system is at the origin. They define 
the scale of the photocentric orbit whose elements are 
those of Baize’s visual orbit except for scale and a rotation 
of 180° in the plane of the sky. The central portion of the 
figure has been rotated 180° with respect to the outer 
portion to permit a direct comparison of the visual and 
photographic observations and to illustrate graphically the 
comparative sizes, a/a = B — 8, of the similar photocentric 
and relative orbits. 


Other trigonometric parallax determinations 
adjusted for the Yale precepts (Jenkins 1952) are 


Allegheny +069 + "008 (p.e.) 


McCormick +.074 + .oII (p.e.). 
-The weighted mean value combined with the 


current Sproul value (unadjusted for Yale pre- 
cepts) yields 
Tabs — 


+"074 = 7003 (p.e.). 


The value of a from the combined solution 
leads to 


QleR 


= B — 6 = +.278 + .007 (p.e.) 
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assuming no error in the value for a. The value 
for Am = 1.83 given by Eggen (1956) leads to | 
B = .156. , 

Adopting the above P, @ and maps, we find: 


a= 8.22 .3.a.u. 
Ma + Mg = 1.12 + .140 


the p.e. arising from the error in a (a.u.) alone. 
The individual masses are 


Ma = .630, Mz = .490 


de Kamp’s compilation of masses (1954). The 
adopted parallax yields absolute visual magni- 
tudes 3.54 and 5.37 for the individual compo-_ 
nents. Using the bolometric corrections adopted _ 
by Eggen (1956), we find absolute bolometric 
magnitudes 3.49 and 5.31 for the A and B com-~ 
ponents respectively. Both stars clearly lie above 


’ the main sequence mass-luminosity curve. This 


departure from expected positions appears greater | 
than their positions on the H-R diagram or color-_ 5 
magnitude diagram compared to main sequence 
stars. 

Permissible changes in parallax, semi-axis ma- 
jor and period can hardly explain this deviation. 
Eggen’s photometric parallax (1956) of "062 
would bring the masses in agreement with the 
mass-luminosity relation of main sequence stars; — 
however, this seems an unlikely value in view of i 
the three trigonometric parallax determinations 
and the small p.e. of their mean. The period — 
cannot be changed sensibly, while an increase of 4 
the order +715 in a, would imply abnormally 
large systematic errors in the visual observations. 
The major axis lies close to the line of nodes, 
which strengthens the value of a derived by — 
Baize. The undermassive quality of both com- 
ponents remains unaccounted for. Further obser- — 
vations of the parallax and orbital motion of this” 
close pair are desired. 1 


=iikineepnmeahb ademas aia 
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PARALLAX AND MASS-RATIO OF FURUHJELM 46 
FROM PHOTOGRAPHS TAKEN WITH THE 24-INCH SPROUL REFRACTOR 


By SARAH LEE LIPPINCOTT 
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q Absiract. Measurement and reduction of photographic plates taken with the Sproul refractor over the interval 1938- 
1958 yield +%1582 + “ood (p.e.) for the relative parallax and +0537 + "0039 for the semi-axis major of the photo- 
lcentric orbit. The adopted value for the absolute parallax, +”%160 + "003, leads to masses .26© and .25© for the primary 


and secondary components respectively. 


Introduction. Furuhjelm 46 = BD +45°2505 
|= HD 155876, 17°9™2, +45°50’ (1900) is a visual 
binary with period 13.02 years (Van Biesbroeck 
1959) and semi-axis major "71 (Baize 1952). The 
‘combined magnitude is 9.38, the magnitude dif- 
ference Am = 0.37 and spectrum of the brighter 
‘component is dM4 (Eggen 1956). The radial 
velocity of the system is —18 km/sec (Wilson 
1953). Recent adjustments of the orbital ele- 
‘ments combined with a 50 per cent increase in 
photographic material justify a new determina- 
tion of the parallax and mass ratio. The current 
‘investigation includes the material used earlier 
and supersedes the former results (Binnendijk 
1951). 

_ Material and measures. The procedure described 
in previous Sproul papers has been followed (van 
‘de Kamp and Lippincott 1949). The material 
includes 860 exposures on 230 plates, taken on 65 
‘nights, resulting in a total weight of 171. The 
‘material is summarized in Table I. The regular 
‘exposure time for one image was about 7 minutes 
\in 1938 and 1939, 14 minutes between 1940-1945 
and 14 minute for 1946 on. The dependences and 
their yearly changes are given in Table II which 
also includes the standard frame and position of 
Furuhjelm 46 for 1940.0. The spectra for the 
reference stars were kindly furnished by Dr. 
Vyssotsky. The image of Furuhjelm 46 was re- 
duced to magnitude 10.0 by the use of a rotating 
‘sector of 58 per cent opening, thus providing 
close magnitude compensation with the three 
reference stars. The plates from 1938 to 1950 
were measured by the earlier investigator (Bin- 
nendijk 1951). The 1951-1958 plates were meas- 
ured and reduced by Mrs. Edward Allen who 
also remeasured the plates taken between 1941 
April 8 and 1943 July 19. The average agree- 
ment between the two measurements is well 
within a micron; the average for each night was 
adopted for further use without increase in weight. 
All plates were measured on the Gaertner Ma- 
chine. 


Solution for parallax and photocentric semi-axts 
major. The maximum separation of the compo- 
nents is I7I or .06 mm; the images appear elon- 
gated under good seeing near the times of maxi- 
mum separation 1944-1946 and 1957-1958. How- 
ever, it has been assumed that the measured 
position of the blended image always represents 
the center of light intensity or the photocenter 
of the two components. The night means are 
represented by the following equation of con- 
ditions: 


Sea (Y¥2) = Co + pet + rPa + aQa 
n= Cy + nyt + Ps + aQs, 


where the symbols have their usual meaning, 
t is counted for 1940.000. The values prior to 
1941 .822 were corrected for a color equation 7z, 
amounting to +.0016 mm due to the adjustment 
of the components of the objective on that date; 
the color equation after 1949.05 amounts to 
+ .0004 mm (Lippincott 1957). 

The orbital factors were computed from Baize’s 
orbit determination (Baize 1958) except for 
changes in period and periastron passage recom- 
mended by Van Biesbroeck (1959) in view of 
subsequent observations: 


Ip 0 aiehoy) 1 = +£153°9 
1 — TOS 2/n18 w= Ee 
C73 OF aaa Fis 0 
phy aes 


Separate solutions in right ascension and declina- 
tion yield the following results: 


p.e. Weight 
Cr = +8.59639 
Cy = +1.86911 
Me = +0.015718 = +0"72966 +0%0005 3981 
Hy = —0.087320 = —1.6477 +0.0006 3955 
tT: = +0.008312 = +0.1568 0.0041 63.91 
Ty = +0.008924 = +0.1684 +0.0125 9.37 
dz = +0.002585 = +0.0488 +0.0050 41.92 
ay, = +0.003161 = +0.0597 +0.0062 38.65 
p.e. Ir = +0.00174 = +0.0328 
p.e. Iy = +0.00205 = +0.0387 
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TABLE I. OBSERVING DATA, MEASURED POSITIONS AND RESIDUALS 
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ke i) Stl athe stl ae 


+ 


Q5 
.67 
.62 
Bh ir) 
-37 
-35 


fe Sea 


Se sinh (le| 


Li +++ 


esheets 


+34 
.18 


x. Vs 
unit ,0oooO! mm 
+8 5829 +2 0194 

5797 = 2, 0225 
5941 19380 
5867 9231 
6047 8329 
+8 6039 +1 8317 
6295 7629 
6267 7581 
6225)" 07387 
6429 6739 
~+8 6389 +1 6673 
6571 5863 
6500 5668 
6736 5037 
6722 4987 
+8 6671 +1 4825 
66901 4807 
6643 4720 
6886 4201 
6872 4100 
+8 6839 +1 4056 
6805 3986 
6832 3982 
7052 3395 
7006 3263 
+8 6991 +1 3119 
6978 3095 
6988 3067 
6972 3015 
6961 2994 
+8 7211 +1 2445 
7172 2388 
7354 1565 
7365 1549 
7281 1420 
ON27 Leal O2 
7287 1374 
7520 0698 
7506 0667 
Mone 0649 
+8 7508 +1 0642 
7508 I 0633 
7669 9788 
7655 9801 
7793 8844 
487716 + 8680 
7713 8659 
7795 8659 
7972 8009 
7911 7888 
+8 8133 + 7108 
8067 7035 
8083 7002 
8347 6253 
8267 6206 


} 
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TABLE I (continued) 


H.A, No. of x ¥ Ux Vy 
Epoch Date min, pl. exp. Obs. Pa Pé Qa 05 unit .0ooOr mm unit .0ooor mm sp 
GG en vec omeaeto. Kr 607, =-3730 -+-1.25 -- .07 4-+-88452 + 5373 —10 — 1 3 
o7Omeviaye10) — 8 “4 16 Fr -+-.405 -.850 1.25 08 8445 5350 — 6 Oye 73) 
Wo7meiiiew Om ite 4°16 Er -+-.064 +-.942 1.26 10 8452 5308 +21 @ 3 
.493 June29 — 8 4 15 Wy —.334 +.904 RAY II 8411 5218 +2 —29 3 
Omelette ta 7 Ko = 3517 -.830 1.28 12 8437 Gey aReO) SP Rt 
534) Julyst4  —22 4 13 Fr —.559 +.806 -+1.28 + .12 +8 8412 + 5195 +16 — 8 2 
eo. 222me vara?) — 0 4 16 )6«6Po 6.966. 186 1.36 -29 8604 4521 =A) sey 
Mo7-405 8 Maye3t —4r 4 16 Wy -+.140 +.932 1.42 56 8779 AOS PAA, eR B 
Se peyetae— 1254 16) Kr —.566 +'.802 1.42 58 8735 BAO La? ae 2 Sane 
Ron{Osmemicnyi2Onse— TO") e4,. 16° Po =-.210  +4-.917 os ap sos 8011 27 Ooo ea 
Bi = Leendert Binnendijk Ke = Geoffrey Keller Ro = Karl Hans Roth 
Da = John J. Davis, Jr. Ko = Robert H. Koch Sl = Morton L. Slater 
D = Roy W. Delaplaine K = Michael S. Kovalenko T = Armstrong Thomas 
Dt = Daniel F. Detwiler m = Peter A. Morris < = Peter van de Kamp 
fl = Edith Flather P = John H. Pitman Wo = Charles E. Worley 
Fr = Laurence W. Fredrick Po = William Poole, Jr. Wy = Arne A. Wyller 
Y = George B. Yntema 
TABLE II, REFERENCE STARS 
i Dep. 
No. Name Mpy Spectrum Diameter Xs Ve I940 AD/yr 
I BD +46°2271 10.2 2 “108 de 62 441.22 220 — .00035 
2 +45 2506 9.3 Go -158 +16.67 — 67.00 294 8 
3 +45 2507 10.3 Go -099 +27.95 +25.78 - 486 Sy | eae 
) Fu 46 (10.0) M4 116 + 8.60 + 1.86 


A combined solution for 7 and a gives the follow- 
ing results: 


p.e Weight 
mm 
Cz = +8.59609 
Cy = +1.86963 
Mz = +0.015721 = +0%29666 +%00056 4019 
My = —0.087319 = —1.64771 =+.00056 4027 
a = +0.008382 = +0.1582 + .0041 73.38 
a = +0.002839 = +0.0537 + .0039 82.57 
p.e. I = +0.00189 = 40.0356 


Table I gives the residuals from the combined 
solution for each night. 

Discussion. Table III gives the mean epoch 
for the normal places of residuals in x and in y, 
together with the total weights and number of 
nights. Considering the weights of the mean re- 
siduals, no significant trend exists. The wide 
range in separations of the components over the 
period has produced no obvious distortion in the 
photocentric orbit. 

The present determination of the relative par- 
allax +1582 + "0041 supersedes the earlier 
Sproul determination which is included in the 
present investigation. The value remains vir- 
tually unaltered. Although the weight of the 
parallax determination has been increased from 
45 to 73, the probable error has not been dimin- 
ished. The p.e. 1, = - "026 for the earlier study 


was artificially low; the average Sproul value is 
+033 (Lippincott 1957). The p.e. 1, for the 
present material yields +07033, and hence a 
more realistic value for the p.e. of the parallax. 

The usual statistical procedure for determining 
the parallax of the reference system (Lippincott 
1958) has been modified due to the large proper 
motion of "132 for reference star no. 2 (Table II). 
Dr. Vyssotsky’s spectra do not permit the classi- 
fication of luminosity type. In this case, the Go 
probably indicates a dwarf star both from the 
point of view of paucity of Go stars of greater 
luminosity and the high tangential velocity indi- 
cated even for the case of a dwarf star. Assuming 
a GoV, a parallax +"o1 is adopted; the average 
parallax for the other two reference stars is found 
in the usual Sproul manner, resulting in a 
weighted mean parallax +0062 for the reference 
system. The value of the absolute parallax of 
Furuhjelm 46 for the Sproul material 1938-1958 
is therefore +’1644 + "0041. 

Other trigonometric parallax determinations 
adjusted for the Yale precepts (Jenkins 1952) are 


McCormick +147 + "013 (p.e.) 
Yerkes + .136 + .o12 
Mt. Wilson + .149 + .O12. 
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TABLE III. NORMAL PLACES OF RESIDUALS 
Vz 'y 
Epoch =p =n unit .00oor mm 
1938.82 4 4 +22 = 6 
40.99 13 5 +10 ear. 
42.90 12 4 — I — 2 
44.44 13 5 ard =e 
45.42 14 5 =r +14 
46.47 19 7 S04, 0 
48.07 20 7 = 9 6 
49.60 19 7 ip co) 
51.80 18 6 = FE — 8 
53.82 14 5 + 2 amc 
55-43 15 6 ap — 6 
57.28 10 4 Bia ais. 
The weighted mean value combined with the 


current Sproul value (unadjusted for Yale pre- 
cepts) yields 


Tabs = +”160 + "003. 


The value of a from the combined solution leads 
to 


= B — 6B = +.076 + .006 (p.e.) 


the probable error arising from the error in the 
a alone. The value given by Eggen (1956) for 
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Am, .37, leads to 6 = .417. Adopting the above 
P,aand taps, we find 

a= 4.44 + .08 a.u. 

Ma, + Mg = -516 .028© ; 


the p.e. arising from the error in a (a.u.) alone. 
The individual masses are 


Ma = 
Ms = 


.262© 
.2540 


where the probable errors are about +.02© due 
to the probable errors in the parallax, the a and 
the Am. The absolute visual magnitudes of the 
two components are I1.01 and 11.39, respec- 
tively. 
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NOTES AND OBSERVATIONS 


NOTE ON ORBITAL ELEMENTS OF DOUBLE STARS 


We may define as the elements of the orbit of 
a double star any set of seven independent pa- 
rameters which, given the double star’s right 
ascension and declination, fix the orientation of 
the orbit plane in space, the orientation of the 
orbit in its plane, the scale of the orbit and the 
motion of the companion along the orbit as a 
function of the time. 

The conventional elements—period P, epoch 
of periastron passage 7, eccentricity e, semi axis 
major a, inclination 7, (ascending) node 2 and 
longitude of periastron w—satisfy these condi- 
tions. The well-known fact that, in the absence 
of observations of the component along the line 
of sight, it is not possible to distinguish between 
ascending and descending node is not really of 
importance in this respect. 

Other systems of elements, which satisfy the 


conditions mentioned, have been suggested. It 
depends entirely on the purpose for which the 
elements are needed which system will be found 
most convenient. If, for example, one wished to 
study the distribution of the orbit planes, the 
substitution of the right ascension and declina- 
tion of the positive pole of the orbit plane for the 
inclination and node would be advantageous. Or, 
if one were interested in computing dynamical 
parallaxes, the substitution of #; = aP- for the 
semi-axis major, and so on. 

It would appear that for most if not all purely 
computational purposes (determination of the 
orbit, ephemeris, differential corrections, etc.) 
the most convenient system so far suggested has 
as the ‘“‘elements’’: the mean motion n, T, e and 
the Thiele-Innes constants A, B, F and G (van 
den Bos 1926). 


| 
I 


1959 November 
__D. Belorizky (1959) has given a graphical 
method for the determination of the ephemeris 
for visual double stars. 

If I understand the author correctly, he there 
advocates the use of Thiele’s constants a, A, b 
and B (Thiele 1883), which he denotes a, p, 0; 
and q, in preference to the Thiele-Innes constants 
A, B, F and G. This would seem to me to be a 
retrograde step, largely because the Thiele-Innes 
constants are functions of the four geometrical 
elements a, 7, w and Q only, while the Thiele 
constants are functions of these four plus the 
eccentricity e. 

Belorizky states: ‘“These elements are thus not 
well adapted for the calculation of the elements 
of the orbit, and particularly for that of the half- 
major axis. Moreover, one has to turn back to 
@;, b;, p and q for the construction of the apparent 
orbit. These last quantities are more convenient 
for the determination of the real orbit (Belorizky 
1957), and appear quite naturally in that of the 
apparent orbit.” 

For those who prefer geometrical construction 
to computing, I have given (van den Bos 1932, 
1933) a method for constructing the apparent 
orbit directly from the data A, B, Fand Gande, 
which seems to me simple enough. The easiest 
way to “‘construct”’ the apparent orbit from orbi- 
tal elements is undoubtedly to plot points from 
the ephemeris, which one has already had to 
compute in any case. 

As far as I can see, in Belorizky’s construction 
drawing MC parallel to OR would not locate the 
point C, unless the apparent orbit had already 
been drawn. Or, alternatively and simpler, one 
would have to construct the length of MC ina 
similar way to that of OM, by drawing a semi- 
circle on RO as diameter. 

For the conversion of the Thiele-Innes con- 
stants we have: 


A+G = a(1i + cos7) cos (w + Q) 
A — G = a(i — cos?) cos (w — Q) 
B— F = a(1 + cos2) sin (w + 2) 
-B — F = a(t — cos7) sin (w — Q) 


The quadrants of w + 2 and w — Q are unam- 
biguously determined by the Thiele-Innes con- 
stants, as a, I + cosz and I — cos7z cannot be 
negative. The quotient of the first and third 
equation determines » + Q, that of the second 
and fourth w — Q. This still gives two solutions 
for w and Q; we take the one where 2 < 180°. 
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Having found w + 2 and w — Q, we have 
a(1 + cost) = (A + G) sec (w + 2) 
= (B — F) cosec (w + Q) 
a(I — cost) = (A — G) sec (w — Q) 
= (—B — F) cosec (w — QQ). 
Adding these gives 2a, after which dividing either 
by a gives cos 7. 


Thiele’s solution of the same problem, using 
Belorizky’s notation, is: 
y cos G = a1 
ysinG = dbisec gy (sin ¢ = @) 
giving the auxiliary quantities y and G. Then 
cos 2I' cos 2w = cos 2G 
cos 2F sin 2 = —sin 2G cos (¢ — p) 
sin 2r = sin 2G sin (¢ — p) 
giving w and the auxiliary quantity [. Then 


cos 2I' cos (¢ + p — 22) = cos (¢ — p) 
cos 2I' sin (¢ + p — 2Q) = cos 2G sin (q — ) 


giving 2, when remembering that w» and p — Q 
must be in the same quadrant. Finally, a and z by 


—171 COS eG Cos 7 — 7. sine Nh. 


I entirely fail to see that this process is simpler 
than the one above, and, therefore, that the 
Thiele constants are better adapted to the calcu- 
lation of the elements of the true orbit than the 
Thiele-Innes constants. 

Belorizky’s solution (1957) is not directly com- 
parable, as he starts from different data, related 
to the axes of the apparent orbit. It therefore 
appears to be distantly related to another method 
(van den Bos 1927), which is also based on the 
axes of the apparent orbit. 

For the computation of an ephemeris, the com- 
parison is: 

Thiele constants 
x = a, cos p(cos E — e) + bi cosqsin E 
y =a, sin p(cos E — e) + bi singsnE 


Thiele-Innes constants 
x=AX+ FY 
y = BX +GY. 


In the first case one has to find E (say, by 
Astrand’s table), sin Z, cos E, cos E — e, sin p, 
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cos p, sing and cosq; in the second, X and Y 
from the X, Y tables. Then, in the first case, six 
multiplications (a; with cos E — e, b; with sin E 
and the two products with cos p and sin #, cos q 
and sin g, respectively) against four in the second 
case. The first process hardly seems simpler than 
the second. Thiele himself gives as the easiest 
process : 


—32(¢@+9)] 
= (b, sin E + a, cos E — aye) cos$(q — fp) 


psin [0 — $(p + g)] 
= (b,sin E — a,cos E + aye) sin #(q — p) 


pcos [1 


which does not seem preferable to the Thiele- 
Innes constants either. 

In differential correction of the orbit, the great 
convenience of the Thiele-Innes constants in 
computing the equations of condition would, I 


believe, be entirely lost by returning to the Thiele. 


constants, apart from the fact that all seven un- 
knowns would appear in the equations, whereas 
in the other case the equation in x does not con- 
tain AB and AG, that in y, AA and AF. This 
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fact obviously shortens the derivation and solu 
tion of the normal equations. 

The above considerations appear to show tha 
the substitution of the Thiele constants for thi 
Thiele-Innes constants would be a retrograd 
step in every respect. 

On one point I find myself 1 in complete agree 
ment with Belorizky. It is that the fact that 
may safely assume that both Thiele and ina 
were familiar with the Gaussian constants (and 
presumably, with Lagrange’ S papers) in no wa} 
detracts from their merit in having shown, fo 
the first time, how to simplify computations ii 


double star orbit work. 
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NOTICE 


The National Science Foundation announces that the next closing date for receipt of proposals fe 
support of renovation and/or construction of graduate level (doctoral) research laboratories is Mare 
I, 1960. Proposals received prior to that date will be reviewed during late spring and early summei 
Disposition of approved proposals will*be made during late summer 1960. Proposals received afte 
the March 1, 1960 closing date will be reviewed following the next closing date which is expected t 
be Geneeniee I, 1960. 

This program will continue to require at least 50 per cent participation by the institution wit 
funds derived from non-Federal sourcés. Proposals may be submitted for modernization or constru 
tion of research laboratories, including laboratory furnishings but not including apparatus or equij 
ment, in any field of the natural sciences. For the present, this Program is restricted to those depar 
ments which have an on-going program leading to the Ph.D degree. Support of facilities to be use 
primarily for instructional purposes will not be considered. It is suggested that, depending on th 
discipline involved, preliminary inquiry be made to either the Division of Biological and Medic: 
Sciences or the Division of Mathematical, Physical, and Engineering Sciences, National Sciene 
Foundation, Washington 25, D. C. Information concerning the Program and instructions for preparé 
tion of proposals may be obtained upon request. 


